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A computer-aided data acquisition system was developed 
and a microthermocouple probe constructed to obtain thermal 
distributions in turbulent buoyant jets exposed to a cross-— 
flowing ambient fluid. The system performed high speed 
temperature measurements as a microthermocouple probe was 
automatically traversed through a sequence of preprogrammed 
positions under the control of a microcomputer. Operability 
of the apparatus was demonstrated by measuring temperature 
distributions in planes perpendicular to the streamwise axis 
of jets from which contour plots of temperature were gener- 
ated. Using temperature distributions along with velocity 
distributions allow buoyant jet characteristics to be 
computed, including the entrainment rate of ambient fluid, 
jet trajectory, and heat transfer to the ambient. The 
experimental technique is discussed and temperature contour 


plots for a jet at various planes are presented. 
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NOMENCLATURE 
PICrementeal Cross-Sectional Area in the 
Temperature Matrix 
Jet Half-width 
Normalized Jet Half-width 
Specific Heat 
Diameter of the Jet at the Nozzle 
Binary Mass Diffusion Coefficient 
Densiometric Froude Number 
Acceleration of Gravity 


Heat Transfer Rate from the Jet to the Ambient 
Fluid 


Ambient-to-Nozzle Flow Ratio 

Length of the Probe Arm 

Radial Distance from the Center of the Jet 
Length of the Probe : 

Schmidt Number 

Streamwise Coordinate Along the Jet Centerline 
Normalized Jet Temperature 

Ambient Fluid Temperature 

Jet Temperature Within the Temperature Matrix 
Centerline Jet Temperature 

Nozzle Temperature 


Jet Temperature As Measured by the Probe 


Temperature Within the Jet at a Radial Distance 
(r) from its Center 


Centerline Velocity of the Jeuyarteene Nozzic 
Discharge Velocity of themger 


Velocity within the Jet at a Radial Distance "r" 
£FYOM Jes eeemrer 


Jet Velocities Corresponding to Locations within 
the Temperactume Marni 


Normalized Centerline Velocity 


Entrainment Coefficient; Offset Angle of the 
Probe Arm 


Offset Angle of the Probe Mounting Bracket 
Probe Angle of Deflection from Horizontal 


Local Angle of Inclinatvon from for tZones om seme 
Jet Streamwise Axis 


Spreading Ratio 

Kinematic Viscosity 

Density of the Jet Fluid 

Density of the Ambient Fluid 

Density of the Ambient Fluid at the Nozzle Exit 
Density of the Jet at centenmmne 


Centerline Density of the Jet at the Nozzle 


Angle of Inclination of the Data Plane from Horizontal 
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I. INTRODUCTION 


Buoyant jets are very common in nature. We see them in 
the form of exhaust gases emitted from smoke stacks of 
refineries, mills and ships. We see them in the form of 
heaved waste water expelled into the sea from power plants 
and from the main propulsion condensers in steam driven 
ships and submarines. It is no wonder that the fluid 
mechanics and heat transfer characteristics of buoyant jets 
have been of interest to environmental, civil and mechanical 
engineers for decades. To evaluate their ecological impact, 
and of most recent interest, to harness buoyant jets as a 
means of detecting military targets and guiding weapons, it 
is necessary to develop models which aerate one prea let 
their trajectory and decay. 

Most studies to date have dealt with buoyant jets rising 
Beeugiea qulescent ambient fluid; however, in nature most 
problems involve flowing ambient fluids. Relatively little 
experimental “on has been done with buoyant jets in cross- 
flow and, according to Hilder [Ref. 1], the trajectories of 
jets and the entrainment rates of ambient fluid predicted 
from previous work do not agree well with one another. Most 
mathematical models of buoyant jets in crossflow assume 
Gaussian profiles for velocity ag temperature. Nickodem 


[Ref. 2] has shown through experiments that in fact, the 


ale 


Gaussian profiles of velocity are altered by crossflow: 
This leads one to suspect that the same may be true for the 
temperature profiles. 

The objective of this work was to develop a system to 
thermally map a buoyant jet in crossflow. Then, by measur-= 
ing both velocity and temperature distributions, improved 
computations of entrainment, trajectory and heat transfer 
characteristics of jets can be made thereby giving rise to 


more accurate models. 


ile 


iM oUOLANT JETS DISCHARGED TO A CROSSFLOW 


ee ROPERTIENS GF BUOYANT JETS 

A buoyant jet is characterized by a momentum and a 
density differential between the jet and its surrounding 
ambient resulting from a variation in temperature and/or 
Fluid concentrations. Therefore, fluid motion in the jet 
is governed by both inertial and buoyant forces. The non- 
dimensional ratio of these forces, known as the densiometric 
Froude number, provides an important quantitative measure- 


ment of jet characteristics and is shown below. 


7 
O 


gD(p_ = oe 


where US is the jet's discharge velocity, g is the accelera- 
miommer gravity, D asethe discharge diameter of the jet, 
oa is the density of the crossflowing ambient and QO, 1s 
wmemdensity Of the jet fluid at its point of discharge. 

The Gaussian velocity and temperature profiles assumed 
by most models of buoyant jets are very similar. Velocity 


behavior 1S given by: 
ULE). a Ue. eee 


where Os to wtheorcenteuiline velocity, © is the independent 


variable and a radial distance from the centerline of the 


ues 


jet, and B is defined as nominal jet halfwidth. As r 
approaches B, velocity decays to (1/e) U_ (Ret: S1tsaeesamnar 


larly, temperature behavior is given by: 


T(x) = T exp(-r°/\*B*) 
where Th 1s the centerline temperature, r and B are defined 
the same as above and i, a spreading ratio, is the inverse 
of the turbulent Schmidt number (S)¥e=s is aerinea aoe emre 
ratio of the molecular momentum and mass diffusivities and 
is equal™= to V/Dap where v is the kinematic viscosity and 
Dap 1s the binary mass diffusion coefficient associated with 
substances A and B [Ref. 4]. Although \A varies inversely 
with the Froude number, the change is very slight, and in 
the case where substances A and B are both water, i is 
slightly greater than 1]. Hirst [Rete Si) found weoeva., 
between 1.16 at F = 0 to llll at F = oeiginde ete dee 
effect then, is a more gradual temperature decay than was 
found with velocity. 

Most buoyant jet models consider the entrainment of the 
ambient fluid into the jet and are based on relevant conser- 
vation equations of mass, momentum and energy. In conserva- 
tion of mass, the downstream change in total mass of the jet 
is equated to the mass of the entrained fluid. The conser- 
vation of momentum must consider both vertical and horizontal 
contributions. Changes in vertical momentum are equated 


to the buoyant forces while changes in the horizontal 
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momentum of the jet are equated to the horizontal momentum 
of the entrained fluid. The conservation of energy involves 
energy changes resulting from variations in the ambient 
temperature as caused by the jet. Hilder [Ref. 1] developed 
the following governing equations in non-dimensional differ- 


ential form. 
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CONTINUITY Ae ee aes 2ab[|u -R cos 6| + a,Rsin 6] 
HORIZONTAL MOMENTUM S=(uAb* cos@) = 4Rob[ |u_—Rcos 6 | +a Rsin0] 
oe Ze 
VERTICAL MOMENTUM aes cee Bee 
ds‘ -m 6 2 
aor mo iD 
he At 
ENERGY sz(u_ T-A 2 = 2 {u_b*} 
(A +1) 


Pee te sOWw. REGIMES 

The jet passes through several regimes as it travels 
From the nozzle through the ambient. The three regions most 
frequently referred to are shown in Figure l. They are the 
zone of flow establishment, the zone of established flow 
and the far-field zone [Ref. 3]. In the zone of flow estab- 
lishment, the velocity and turbuleygce profiles transform from 
the conditions within the nozzle to a free turbulent flow 
condition. It is in this region that the jet begins to mix 


with the ambient fluid; however, the flow is still more 
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pemougliye ine luenced by the nozzle discharge conditions than 
by the ambient. When the turbulent mixing has reached the 
Pemeerline of the yet, the zone of established flow is said 
eo begin.» in this region, the profiles have assumed their 
free turbulent shapes. Now the flow is governed by the 
jets' momentum and buoyancy as well as by the condition of 
the crossflow. The far field zone is defined as that region 
in which jet momentum is depleted and the jet fluid is con- 


vected and diffused by the ambient currents and turbulence. 


Cees riers OF CROSSFLOW 

AGetme iammediate exit Of a cylindrical nozzle, a verti- 
cally discharged buoyant jet has a nearly uniform velocity 
distribution and has the same cross-sectional shape as the 
nozzle itself. The velocity gradient between the jet and 
the crossflowing ambient creates longitudinal shear stresses 
at the jet's sides, a positive pressure region immediately 
upstream and a negative pressure region immediately downstream 
of the jet. This results in the deflection of the jet's 
trajectory in the aernserean Grreeceions(Pugure 2), the 
Saeattom OF COUuNnterrotating vortices at the jet's outer 
edges and the deformation of the original circular cross- 
SBecetonal shape into the form of a kidney. As the stream- 
wise axis of the jet approaches the direction of the cross- 


flow, these effects become progressively less pronounced. 
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Ill. EXPERIMENTAL APPARATUS 


eee oYoLEM OVERVIEW 

A surplus milling machine was configured with synchron- 
Popo ive merors amterfaced with a microcomputer that auto- 
matically positioned its bed. It was used as a three-dimensional 
positioning platform in the same manner as in the laser 
Doppler velocimetry work undertaken by Nickodem [Ref. 2]. 
The milling machine was placed adjacent to a rectangular 
plexiglass flume through which the crossflowing ambient 
fluid flowed. A vertical nozzle was installed in the base 
Demene feune eo provide the jet. A temperature probe was 
suspended through an opening in the top of the flume above 
the nozzle by an arm attached to a base mounted on the 
milling machine bed as shown in Figure 3. As the probe was 
automatically traversed through a series of preprogrammed 
positions across the jet, temperature data was automatically 


sensed and stored at high speeds by the computer. 


Bee ChOgo- LOW SYSTEM 

As illustrated in Figure 4, the crossflow circulation 
Pump LOGI water from the cylindrical 248.8 1 (65.7 gal) 
reservoir shown in Figure 5 and discharged through 5.076 cm 
(2 in) diameter tubing into a cylindrical flow settling 
chamber 30.46 cm (12 in) in diameter and 60.91 cm (24 in) 


Peeleeatedwwathim the 60.91 cm x60.91 cm x 88.83 cm 
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eee 251 x35 in) inlet chamber shown in Figure 6. The 
settling chamber was sealed at its bottom so that the water 
spilled from its top into the inlet chamber through honey- 
combed flow straighteners to reduce turbulence and evenly 
disperse the flow. To further reduce turbulence, the flow 
was broken by another stack of honeycombed flow straighteners 
and a layer of fiberglass filter material located immediately 
above the normal operating water level. The flow next 

Ce eee csmcm ~~ o2.59 em X1l3e2.9 cm (9.625 inx12.75 inx72 in) 
flume shown in Figure 7 through a vertical section of the 
same honeycombed material mentioned above. To avoid inadver- 
tent spillage over the sides of the flume during system 
start-up, a 5.076 cm (2 in) diameter overflow pipe was 
located in the inlet chamber. During normal operation, 

a gate valve in this piping was closed. The flow left oe 
flume through a 7.614 cm (3 in) diameter pipe at its end 

and re-entered the crossflow circulation pump reservoir. A 
gate valve located in this piping and shown in Figure 8 was 
used to regulate the water level and flow velocity in the 
flume. The optimum adjustment of this valve was determined 
by trial and error to be closed two turns from its fully 
Open position. Either a globe or ball valve would have been 
more appropriate for this purpose; however, neither was 
readily able, so the gate valve was used. The bracket 
shown at the base of the flume in Figure 8 maintained align- 


ment between the flume and the milling machine. The water 
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Figure 7. 


Figure 8. 


Flume Arrangement 
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Flume Discharge Piping 
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in the flume was gradually heated by repetitive circulation 
through the crossflow pump and by the addition of the heated 
water from the jet. To maintain a constant temperature 
eressrtlow, cooling mater from a refrigerated bath shown in 
Reiieeme was Circulated through a coil of 1.269 cm (.5 in) 
diameter copper tubing located in the flow settling chamber. 
Also, fresh water was added at the flow settling chamber 

as an equal amount was drained from the crossflow pump 
BeserVOmomeanougneas).2609 cm (.5 in) diameter pipe. Cross- 
flow temperature was monitored by a Type-T thermocouple 


Located aneene inlet chamber. 


Cc. JET yoru 

In reference to Figures 4 and 5, the jet flow pump 
C1 ECU lave eo bmmmomeaerectangulars26.2Z2771 (6.94 gal) 
reservoir and discharged through 1.26 cm (.5 in) diameter 
poo end  eem—r Ome nmm 0 omen < 40nNem (13 an x20 in x 21.5 in) 
head tank (Figure 10). The amount of flow to the head tank 
was regulated by a globe valve. Due to a low flow rate 
to the head tank, water was also recirculated back to the 
Be=servyorrm If Orders tommalntain SULLICiene f£row through the 
jet pump to prevent overheating it. A constant water level 
waS maintained in the head tank by a stand pipe which allowed 
overflow back to the reservoir. Sufficient flow into the 
tank waS maintained to make sure that it slightly overflowed 
EOvermueusliy., Water drained from the bottom of the head 


tank through 1.26 cm (.5 in) diameter tubing and passed 
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Figure 10. Head Tank 


ou 


through a rotometer, a .95 cm (.3/75 in) tubing reducer, a 
dye injection system (Figure ll), a water heater (Figure 12) 
consisting of approximately 6.09 m1 )(Z0Mi esc. Seer = 
in) diameter copper tubing coiled ina heated bath and 
finally a 7.144 cm (.28125 in) nozzle which discharged into 
the bottom of the flume. Flow was controlled by pinching 
the tubing between the heater and the nozzle with surgical 
clamps. Drainage from the crossflow reservoir discussed in 
Section III.B was used to replenish the jet reservoir. The 
dye injection system, used in photographing the jet, was 
located approximately 8.23 m (27 ft) upstream of the nozzle 
to minimize any disturbance to the jet that it might have 
caused. The majority of this distance was taken up by the 
heating coil mentioned above. The vertical distance between 
the top of the stand pipe in the head tank and the tip of 
the nozzle in the flume was 2.2 m (86.5 in) which equated to 
21.56 KPa (3.127 psig). Jet flow temperature was monitored 
by a Type-T thermocouple located within the jet flow tubing 


approximately 1.167 m (46 in) from the nozzle. 


D. TEMPERATURE PROBE 

Measuring temperatures in a buoyant jet with a thermo- 
couple is intrusive. To reduce the probability of distorting 
results, steps were taken to minimize the cross-sectional 
area of the temperature measuring device as seen by the flow 
of the jet. A .0254 mm (.001 in) diameter Type-E micro- 


thermocouple was selected. The suspension device for the 
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Jet Heater 
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Figure 


Dye Injection System 
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Figure 


microthermocouple had to be rigid and have a small cross- 
sectional area, for reasons discussed above, as well as be 
an electrical insulator to prevent interference with the 
thermocouple performance. A glass annulus approximately 
1.45 mm (.057 in) in diameter and 11.27 cm (4.4375 in) in 
length was chosen. One lead of the thermocouple was 
threaded through the annulus and the other was glued with a 
fast drying modelers’ glue along the outer surface, allowing 
the microthermocouple junction to protrude slightly from the 
tip of the annulus. The leads at the opposite end of the 
annulus were welded to .0762 mm (.003 in) diameter wire 


which subsequently was connected to 28 AWG extension wire 


to the computer. The annulus was mounted as shown in Figure 
13. Henceforth, this device will be referred to as the 
probe. 


E. PROBE ACTUATOR ASSEMBLY 

The cross sectional area of the probe as seen by the jet 
was further reduced by orienting the probe tangentially to 
the trajectory of the jet as shown in Figure 14. This 
photograph indicated that the probe created no noticeable 
interference with the jet hydrodynamics. Probe orientation 
was accomplished by the linkage assembly shown in Figure l5. 
The fixed end of the probe was hinged to a streamlined tube 
23.495 cm (9.25 in) long with a maximum wider anoscc peas 
seen by the jet, of 3.175 mm (.125 in) Sandee. 2 > me oe) 


respectively. It was rigidly connected to the mounting 
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Figure 14. Probe in Jet 
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Drive Motor Pulley 
Cable and Turnbuckle 
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Linkage Rod 
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Streamlined Tube 
Hinge Assembly 

Probe (Glass Annulus) 





Figure 15. Probe Assembly Linkage 
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bracket as shown in Figure 15 which was connected to the 
probe arm shown in Figure 3 by a single stud which allowed 
pivoting of the probe from side-to-side. The stud was also 
hollowed so that a linkage rod could extend from the hinge 
assembly through the tube and stud to a stroking pulley 
which was rotated by a small motor. The hinge assembly and 
the stroking pulley were spring loaded to reduce hysteresis. 
As shown in Figures 16 and 17, the 1.5 VDC motor, geared to 
one rpm, was directly coupled to a potentiometer as well as 
the drive pulley. The potentiometer was configured in a 
voltage divider such that the amount of motor rotation, and 
ultimately the degree of probe deflection, was proportional 
to the potential difference sensed across the potentiometer. 
Limit switches were installed at the stroke pulley as shown 
in Figures 16 and 18 to prevent damage to the linkage 


assembly due to over-rotation. 


F. MICROCOMPUTER INTERFACE 

The data collection process consisted of adjusting the 
probe angle of deflection, traversing the three-dimensional 
positioning platform and measuring temperature profiles. 
All of the mechanisms which controlled these events were 
interfaced to an HP-9826 computer shown in Figure 19 through 
an HP-6942A multiprogrammer which performed high speed analog- 
to-digital conversions and ultimately provided control 
Signals to govern relays within the system. Refer to MAIN T 
in Appendix B for the microcomputer software which directed 


this process. 
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Figure 17. Probe Actuator Motor- 
Potentiometer Arrangement 





Figure 18. Probe Actuator 
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ee Sooeeenele Adjustment 
As discussed in Section III.E, a potentiometer con- 
figured as a voltage divider provided probe angle feedback 
to the computer as shown in Figure 20. The direction of 
motor rotation was controlled by the computer through Type- 
MRRICDL replays connected as shown in Figure 21. When the 
probe was at a desired position, a 5.0 VDC signal was applied 
to pins 6L and 6R which allowed both relays to assume the 
normally closed (NC) positions which opened the power circuit 
to the motor. When it was desired to rotate the motor clock- 
wise, pin 6R was grounded which resulted in 5.0 VDC applied 
across the coil in the right-hand relay. This caused the 
relay to assume its normally open (NO) position resulting in 
a 1.5 VDC signal at terminal B of the motor cauSing it to 
rotate in the clockwise direction. The left-hand relay was 
activated in a similar manner for counterclockwise rotation. 
Pee tol ROBE SUBS Im Appendix B for probe positioning 
software. 
2. 3-D Positioning Platform Movement 
Positioning platform movement was controlled ina 
manner similar to the probe and was discussed in detail by 
Nickodem [Ref. 2]. Refer to MTR SUBS in Appendix B for 
associated HP-9826 software. 
3. Temperature Data Collection 
Three thermocouples were monitored in the data col- 
iiceronearocess. A Type-T thermocouple located in the inlet 


chamber measured the ambient fluid temperature in the flume, 
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a Type-T thermcouple located in the tubing between the heater 
and the nozzle measured nozzle temperature and a Type-E 
microthermocouple in the probe measured the temperature in 
the jet. The EMF's generated by these thermocouples were 
amplified by "Omega Omni-Amp IIB" millivolt amplifiers shown 
in: Figure, 22 pricrase entering the multiprogrammer for 
analog-to-digital conversion and eventual transformation to 
temperature readings. Fourth-order least squares coeffi- 
cients for this conversion were taken from Beckwith [Ref. 5]. 
Operation of the crossflow circulation pump created suffi- 
cient electrical interference to distort the thermocouple 
Signals. This problem was corrected by applying a thin 
coating of silicon sealant to the Type-T thermocouple junc- 
tions and by connecting the crossflow circulation pump 
casing, the nozzle and the jet tubing in the vicinity of the 
nozzle thermocouple to a common ground. Because the jet 
tubing was plastic, it was necessary to manufacture a brass 
"T" connector as shown in Figure 23 which was grounded and 
located in close proximity to the thermocouple junction. 


Refer to T SUBS in Appendix 8) tem tas seca dso ely cause 
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Figure 22. Thermocouple Amplifiers and HP-6942A 
Multiprogrammer 
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Figure 23. Nozzle Thermocouple Grounding Arrangement 
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IV. EXPERIMENTAL PROCEDURES 


A. CALIBRATION 

Eerore —~nme data collection process could begin, the 
rotometer, the thermocouples, the probe and the positioning 
platform had to be calibrated. The detailed steps taken 
are discussed below. 

l. Rotometer 

With a constant level maintained in the head tank 

and the jet tubing disconnected from the nozzle and elevated 
to the same height as the top of the nozzle, five 100 ml 
samples were drawn through the rotometer and timed to the 
nearest 0.01 second at each rotometer reading from 102 
to 75% in 5% increments. Flowrates and standard deviations 
in ml/s are shown in Table l. 

2. Thermocouples 

Since nozzle and probe temperatures were to be 

normalized by the ambient temperature, the nozzle and probe 
thermocouples were calibrated relative to the ambient thermo- 
couple by using the microcomputer program T CAL in Appendix 
B. The procedure followed is outlined in the initial com- 
ments of the program. Coefficients for first order curve 
fits were solved by the least squares method with the 


mainframe programs TCAL and TFIT found in Appendix C. 
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3. Probe 
As the probe assembly was being developed, it was 
convenient to test its suitability with the probe calibration 
panel shown in Figure 24. Resistance changes across the 
potentiometer were recorded for varying degrees of deflec- 
tion. Analysis of this information led to improved designs 
from the standpoint of reduced hysteresis and repeatability. 
The microcomputer program PROBE CAL in Appendix B was 
developed to enable calibration of the final design after 
it was installed in the system as shown in Figures 25 and 
26. The calibration procedure is outlined in the preliminary 
comments of the program. 
4. 3=-D Positioning Plaeseum 
The positioning platform was calibrated in a manner 
that placed the tip of the probe at desired locations within 
the flume relative to the tip of the nozzle. Referring to 
the coordinate system illustrated in Figure 2, the center 
of the nozzle was defined as (0,0,0) in xyz-coordinates. 
The following relationships apply to the probe geometry 


shown in Figure 27: 


X(real) = X, 7 R, cosa tan(1/4 - a/2) 
Y(real) = - + aoe =—_COS yaa R, cos 4 
Z(real) = Z6 = a Sin y 
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Figure 26. Probe Calibration Panel 
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Figure 27. Probe Positioning Geometry 
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(X01 0120) was the position of the milling machine bed 

when the tip of the probe was at position (0,0,0) with 

y = 0 and a = 8 = 90 degrees. For calibration, the probe 
and probe arm were configured with these settings as shown 
in Figure 28. With measured values cE ro and R, entered 
into the microcompute® program MAIN T7 the weal bia tren a 
accomplished by the program MOTOR CAl theo ec)— yee 
procedure followed was outlined in the subprogram "SUB 
Calibrate." As illustrated in Figures 3, 27 and 28, the 
beng thmens R, could be modified to compensate for adjustments 
of a and 8 to positions other than 90 degrees. Decreasing 
a@ increased the distance along the Y-axis in which the 
probe could be positioned. Increments of a and 8 were 
scribed on the top of the base and at the tip of the probe 
arm in Figure 28 to accommodate this change, if desired. 
The program MAIN T queried Sehne Weems eOmee mc ensee emma cic 
assumed a= 8. The calibration software also established 
position limits to prevent driving the probe into the sides 


of the flume. 


B. PRELIMENAR wes 

Crossflow velocity was determined by injecting blue food 
coloring into the flow and €iming 1es@epavel tnreugheae slam 
interval. The average of several trials indicated the 
velocity was .130 m/s (.427 ft/sec) with the Elume oucter 


valve closed two turns from its fully open position. 
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Figure 29. Typical Buoyant Jet as Obsered with Dye 
Injected 
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Photographs of jet profiles as shown in Figure 29 were 
taken to determine jet trajectory and halfwidths along the 
streamwise axis. This was done by injecting blue food 
coloring into the jet flow as #a@wsceussed in S-eerioneamc- 

The specific gravity of the food coloring was found to be 
considerably less than that of water. To eliminate the 
added buoyant effect this would have had on the jet, a small 
quantity of alcohol was mixed with the food coloring as 
suggested by Merzkirch [Ref. 6]. The amount of alcohol added 
was determined by trial and error. As small quantities were 
added and mixed, samples were gently placed on the surface 
of a beaker of water. Pure food coloring laid on ehereun— 
face and very slowly mixed with the water. As alcohol was 
added, this buoyant effect grew progressively less and the 
mixture would settle into the water. The mixture was con- 
Sidered satisfactory when it no longer laid on the surface, 
but settled to some equilibrium position in the beaker. 

Slide photographs of the jets were projected onto large 
sheets of 3.175 mm (.125 in) grid graph papereand argmri zed 
along approximate streamwise axes and half-width trajectories. 
A scaling factor was determined by equating the projected 
width of the nozzle to its known outer diameter of 7.9375 mm 
(.3125 in). The above data was fit to the Michaelis-—Menter 
Equation [Ref. 7] shown below by the least squares method 
with the mainframe program JETCURV in Appendix C: 


Zi = _ ay 
bate, 
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SCouretetten coefficients elose to 1 were consistently 
obtained. To determine positions within the jet at which 
to make temperature measurements, five evenly spaced posi- 
tions per jet flow rate were selected along the streamwise 
axis in the zone of established flow. Data planes slightly 
larger than the jet width were centered at these points and 
Oriented perpendicular to the streamwise axis. One hundred 
data points were selected in a symmetric square matrix with 
points most densely populated near the center. The planes 
were identified alphabetically and in consecutive order from 
"A" to "F", where "A" represented the plane nearest the 
nozzle. The positions were entered into the microcomputer 
“Mamoronee Dy eptane on a Lloppy disk by the program LOAD XYZ 
in Appendix B. Accompanying each data point was a probe 
deflection angle used to orient the probe parallel to the 
path of the jet to minimize interference. This angle was 
determined by evaluating the first derivatives of the 
equations developed for the streamwise axis and the half- 
widths and performing an interpolation based upon the data 


points' position relative to the two curves. 


Ge DATA ACOUISITION 

The flow systems were placed into operation and the 
ambient and nozzle temperatures were monitored with the 
program T SUBS to evaluate system stability and readiness 
for data acquisition. The system usually took approximately 


two hours to come into equilibrium. This could be monitored 


ak 


by watching the jet nozzle temperatunge in enecen i emen. 
were stable, data consisting of two hundred probe, ten 
ambient and ten nozzle temperature samples per position was 
collected, one plane at a time? by Sone @eregr ane) a ee 

The following information was stored on a floppy disk for 
each data point: “xy. emda :z coordinates; mean probe, 
nozzle and ambient temperatures and the standard deviation 
of the probe measurements. The data was transferred to 

the mainframe computer by using a modem, the microcomputer 


program SEND DATA and the mainframe program GRAB. 


D. DATA REDUCTION 

The raw data was organized into a more usable format 
by the mainframe program TDATA which also converted the XYZ 
coordinates into the XSW system shown in Figure 2. The 
resulting data was selectively sent to the program CONTOUR4 
which applied calibration coefficients to the temperature 


data and normalized it in the following manner: 


where a was the jet temperature as measured by the probe, 
Ls was the ambient fluid temperature and i was the 
temperature of the jet within the nozzle. 

Contour plots of this information, generated by the 
CONTOUR option of the graphics package DISSPLA [Ref. 8], 


are presented in Figures 30 through 35. 


oe 
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Figure 30. Plane A Temperature Contour Plot 
(large scale) 
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Figure 31. Plane A Temperature Contour Plot 
(small scale) 
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Figure 32. Plane B Temperature Contour Plot 


55 








TEMPERATURE GONTOURS IN A BUGYANT JET 
With A GkOSSFLOWING ALIBIED 


30 


20 26 


a 
-10 0 16 


20 





26 


-30 
= 


-30 -26 -20 -16 -10 -6 


Figure 33. Plane C Temperature COmEOUE EE Tee 
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Figure 34. Plane D Temperature Contour Plot 
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Figure 35. Plane E Temperature ContoumePiot 
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fee ULTS 
iemcontour plots found in Figures 30 through 35 support 
PMemeomwyectures in Section I1.C concerning the effects of 
the crossflow. Figure 30, a plot of plane A centered on 
the streamwise axis 46 degrees from horizontal and 7.327 mm 
(.288 in) downstream of the nozzle shows significant distor- 
tion. Figures 31 through 35’ show planes A through E sequen- 
tially plotted on the same scale in order to observe overall 
jet behavior. Plane E was located 87.313 mm (3.438 in) 
downstream of the nozzle and 86 degrees from horizontal. 
it can be seen that as the jet traveled further downstream, 
the distorting effect grew progressively less, as expected. 
The rate of heat transfer from the jet to the ambient 
was calculated for each plane utilizing the temperature 
distribution matrix generated in the program CONTOUR4 in 


Appendix C and the following relationship: 


™m ry 


D = EEE mci 
Be 2 yey? Sp May, MH TA 


where o was the relative density of the jet, iG was the area 
of each matrix segment, Us was the velocity in each segment, 
es was the specific heat of the jet and ae was the tempera- 
ture in each segment. Velocity was measured along the 
streamwise axis by a laser Doppler velocimeter (LDV). It 
pippeated to be nearly constant in the region of the jet 


observed. The mean and standard deviation was 44.875 mm/s 
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and 1.8 mm/s respectively. By assuming the Gaussian profile 
shown in Section II.B, velocity was determined at radii 
corresponding to each segment in the matrix mentioned 


above. The rates of heat transfer are shown in Table 2. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The objective of this thesis was to develop a computer- 
aided data acquisition system and construct a microthermo- 
eouple probe to be used by follow-on students. to study 
temperature distributions in turbulent buoyant jets. Sample 
data was taken to verify system operability. Based on 
results, the system performed in a satisfactory manner and 
will be an invaluable tool for subsequent studies of buoyant 
jets in a crossflowing ambient. 

Data point positions were hand calculated and loaded 
PieOnwtie MiLCrOComputer with the program LOAD XYZ. This was 
an extremely time consuming task and distracted the user 
from defining more than 100 data points per plane. The 
system can be greatly improved with the addition of micro- 
computer software that would automatically determine and 
Weaeedata pOlnt positions. Data point population could then 
be increased with ease which should result in smoother con- 
tour plots and more accurate heat transfer calculations. 

ft was necessary to continuously add fresh water to the 
crossflow system as an equal amount was drained from it in 
order to maintain the crossflowing ambient at a constant 
temperature. This was because the cooling coil located 
within the flow settling chamber was inadequate to compen- 


maeewcOor the heat added to the system by the jet and the 


Gel: 


crossflow circulation pump by itself. Although it was 
possible to maintain the ambient constant within 1.4 € by 
this method, in the interest of conserving water, it is 
recommended that either a larger capacity chilled water 
bath or a cooling system that circulates a refrigerant 


rather than water be appropriated for this purpose. 
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APPENDIX A 


UNCERTAINTY ANALYSIS 


Experimental uncertainty was analyzed in accordance with 
the guidelines set forth by Holman [Ref. 9]. Uncertainties 
in the primary measurements, based upon manufacturer speci- 
Peeditiens and/Or the number of significant digits which 


could be read, follow: 


1. Time (rotometer calibration): .Ol s 

2. Volume (rotometer calibration): iO tole 

3. Rotometer reading: Oo percent 
4, 3-D positioning platform resolution: .1524 mm 

5. Probe deflection angle resolution: .5 degrees 
6. Thermocouple resolution: BOO Oe. C 

Pw elonecrmocouplLe time constant: <0042s5 

fee Lime: (crosstlow velocity) : Poles 

9. Length (crossflow velocity): 705910 


Based on the above, the uncertainty of the flow rate of the 
jet was estimated to be 1.0 ml/s. Maximum uncertainty in 
the position of the probe's tip due to the uncertainty in 
deflection angle was determined to be .983 mm (.0387 in). 
Combined with the resolution of the positioning platform, 
the tip of the probe was positioned with an uncertainty of 
1.135 mm (.0447 in) in each plane, or with an overall 


uncertainty of 1.605 mm (.063 in) in three-dimensional 
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space. Temperature was measured at the approximate rate of 
100 samples per second, well within the constraints of the 

thermocouple time .constant. Uncertainty in the temperature 
measurements were governed by the resolution of the analog- 
to-digital converter which was .0099 °C. Uncertainty in 


the crossflow velocity was .00625 m/s (.0205 ft/sec). 
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APPENDIX B 


MICROCOMPUTER PROGRAMS 


are MAIN T 
10 ' MAIN_T 
9 ? 
3 ' This program coordinates the entire 
40 ! data-taking evolution for measuring 
20 : temperature distributions in buoyant 
60 ! jets. 
70 ! 
80 ! 
30 bee 1. Load ail subprograms 
100 ! 2. Input desired positions from a 
110 i disk file of the form : “RUNEX” 
120 ! (value "XX" is the run number) 
sd ‘ 3. Move the 3-D positioning 
140 ' platform to each position 
150 ! 4. Align the thermocouple probe with 
160 ! the jet streamuise axis 
170 ' 9. Obtain 200 temperatures at each 
180 ! position and compute the mean and 
190 ! standard deviation (sd) 
200 ! 6. Write to disk: 
210 ! AteeCC) 
220 ! b. sd 
230 ! Cc. Xsy¥ ee ein (mm) 
240 ! nozzle centerline is (0,0.0) 
250 ' 
260 ! Probe and arm dimensions in inches: 
2/0 ! aawongntest arm lengtnmame) .0 
280 ! 6b. Longest arm length = 24.0 
290 ' 
300 Length_probe=4, 4375 
310 Length_arm=20.0 
320 ! 
330 OPTION BASE 1 
340 Ditecoaet (12044 (500), 1 Couuy., 20500) ,Probe_angle(500> 
350 LOADSUB ALL FROM °"'T_ SUBS” 
360 LOADSUB ALL FROM “MTR_SUBS"” 
370 LOADSUB ALL FROM “PROBE_SUBS” 
380 CALL Retrieve_coef  (Coef(*),**motor_coef’*’’) 
S50) : 
ny 2. Input desired positions from disk 
0 ' 
420 =| ; 
430 INPUT “Angle of arm relative to +Y-axis?”,Angle_arm 
440 INPUT "Filename for positioning data?"”,Filename’ 
450 ASSIGN @File4 TQ FilenameS 
460 Go_on: ! 
470 ENTER @File4:X¢(I),Y¢I).Z¢1) 
48) - ' 
490 IF X¢€I)<>-100 THEN 
500 IF X¢(1I)2-999 THEN 
a0 Probe_angle(I)=Y¢I) 
520 Prana. em) i> 
oO ELSE 
540 Probe_angle(I)=P_angle 
550 L=1+] 
560 END IF 
570 ' 
580 GOTO Go_on 
530 END IF 
500 ' 
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poate 
620 
630 
649 
530 
660 
670 
680 
690 
700 
710 
720 
730 
740 
70 
760 
770 
780 
790 
800 
810 
820 
830 


t 


GOTO Go_on 
END IF 
' 


ASSIGN oFile4 TO = 
Nitems=I-1 
' 


'3. Begin loop to take data at each point 
' 


INPUT “NAME OF FILE WHERE DATA IS TO BE STORED?” ,FilenamelS 
Records=(Nitems#*8*#/7/256) +2 

CREATE BDAT Filename!t$,Records 

ASSIGN ®@®Filei TO Frilenamel$S 

' 


' 

' 

FOR IT_position=1!1 TO Nitems 
' 





‘a. Move milling machine and move the 


! thermocouple probe. 
' 


1=I_position 
CALL Move_Idv_to(X¢1I) ,Y¥¢1I),2Z¢1),Lenght_arm,Length_probe,Angle_arm, 


e_ angle(I) ,Coef(*)) 
840 ' 


850 
860 
870 
380 
890 
900 
910 
920 
330 
940 
950 
960 
370 
980 
J20 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
Te 
1120 END 


’ 
'b. Obtain mean temperature and sd 
' 


CALL T_couple(T, "PROBE" ,"C",200,St_dev) 
' 


'c. Measure ambient and jet temps 
' 


CALL T_couple(T_ambient, “AMBIENT’,"C",10,Sd) 
CALL T_couple(l_nozzle. NUZ2RE = © tesa) 


'd. Write all information to disk 
' 


QUIPUT@@Fi lelsxcC lo yt) 2 
OUTPUT @Filel;T,T_ambient,T7_nozzle 
QUTPUT @Filel:St_dev 

NEXT I_position 

i 


14, Close files 
' 


OUTPUT @Filel;-100 
ASSIGN @Filel TO « 
f 


1S, Terminate program 
i 

PRINT “All done!" 

B 


Sie 
BEE 
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Zee 


S50 
340 


430 


470 


530 
Sal 


PROBE SUBS 


PROBE SUBS 


This program moves the temperature probe 
to desired angles of deflection. 


! 
' 
! 
! NOTE: Calibration Goeteritc tents are 
! entered in SUB Read_angle. 

' beginning at line 299. 


! 
CALL Read_angle(Angle) 


PRINT) “Ihe erebe is pudsently at ‘Angle; degrees from horizontal. ~ 


INPUT What is your desinied angle for the probe?” ,Desired_angle 
CALL Probe_move(Desired_angle) 
Goi, 120 
END 
' 
! 
' 
SUB Read_angle(Actual_angle) 
(This program reads the present angle 
ToOnetne Probesand Teturns 1t. 
' Q degrees = horizontal 
190 degrees = vertical, downward 
! 
‘AD, Al and A2 = coefficients for a 
!second order curve fit of mV vs angle 
fof deflection data. 
' 


les iatenant oar ommomelunemda taupe lou: 
’ 


A0D=-23.8657824056449 
Al=3.09995506283164E-2 
A= ie eco 2c ooo Ie 6 
' ' FORMAT A/D CARD 
Pewee 2s. CC Ml! Clear Relay Card 
OU 2o 66 .3| 7 Clear ma) Card 
DU coe Coy lo = Oloaripiagi tal: Card 
OUTPUT 7235 25 o.oo come |. 
OWUIPUT £23. 0Bw iO.) eemSE RELAY 
OUMPU Te? Zone ews ! START A/D 
ENTER ee7 2o0iay 
Actual_angle=AD+tA1*V+tA2VeV 
SUBEND 
’ 


! 
! 
SUB Probe_move(Desired_angle) 
'This subprogram moves the probe to the 
‘desired angle 
‘ 
' 0. Check to see if the angle ts in 
! an acceptable range. 
' 


GCUTPUI= (233 "0F 1 Ol |! Glearuke lays 
CALL Clear_screen 
! 


IF Desired_angle>90 THEN 


BEEP 3400.1 
Bier Sa O31 


7 


600 
610 
520 
630 
b40 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
(20 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
ano 
320 
330 
940 
Bo0 
960 
970 
380 
330 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
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1 
1 
1 
| 
1 
| 
1 


— —_ 2 —-) —2) —e 3) —3) — 3 —F 
DW O™N NU & WN 
ooo ooceoc mo 


o——- Om sm Ow ommp om sw 


oD om om oe 


' 

! 

! 
Repeat: 


BEEP 3400,1! 
PRINT “Desired angle exceeds 90 degrees!!!" 
SUBEXIT 
END IF 
' 


IF Desired_angle<0 THEN 
BEEP 3400, 1 
BEEP 3800,1 
BEEP ooU 0g 
PRINT “Desired angle is negative!!!" 
SUBEXIT 
ENOT Tr 
1 Clear the digital] output card. 
- Format the A/D card. 
~ Close the relay that corrects 
the probe potentiometer to the 
A/D converter. 


GUEPUIN? 23 OR 0 

OUTPUT 223: 5R 53585 oy eee le 
GUIBU) 2d Cis 

OUTROS 2s. 05.2 10 


oop 


2. Define the acceptable tolerance 
in the angle (degrees). 


Tolerance=.5 
3. Control woop. 


CALL Read_angle(Actual_angle) 
PRINT “ANGLE =": 
PRINT USING “DDD.DD’:Actual_angle 
BEEP Actual_angle«100,.05 
Angle_error=(Desired_angle-Actual_angle) 


IF ABS(Angle_error)>Tolerance THEN 
' 


ig Angle_error>=0 THEN 
Directions$="Douwn” 
EESE 


DirectionS="Up” 
ENOa Te 
' 


CALL Motor_go(Direction$) 
' 


GOTO penecin 
END IF 


GUIPUT 722 gee ue 
1100 OUTPUT 7232 GC. 11° §'Clearmrelay -eare- 


SUBEND 
' 


i 
i 
SUB Motor_go(Direction$) 
IF DirectionS="Up” THEN 
LUbat=38 


END IF 
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1200 Le eDirect ton>= Down 9 ney 
G20 Port=/ 

1220 ene 

1230 t 

1240 ae «f23% OP 7 22° bpite ol 
1250 SUBEND 

1260 ! 

1270 ! 

1280 ! 

1290 SUB Clear_screen 

Ise aeGrear the CRI. 

Os! 

1320 OUTPUT 2 USING “#.B":255,75 
1330 GCLEAR 

1340 SUBEND 
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3. MTR SUBS 


240 
250 


930 


500 
610 


ATR2SUbs 


4 
' The following series ot subroutines 

f are utilized to calibrate the positioning 
' platform and ultimately to move the probe 
' tip to desired positions within tne jet. 


' 


SUB Draw _f lume 

t £ 

’ Draw the Buoyant Jet Flume on the CRT. 
‘ 


GCLEAR 
GRAPHICS ON 

WINDOW 0,48,0,38 
LINE TYPE 1 

MOVE 6,5 

’ Draw the top view 
IDRAW 
IDRAW 
IDRAW 
IDRAW 
' 


IMOVE 
IDRAW 
IMDVE 
IDRAW 

' lLabe 


Or 
Do 


OnNonN ees Co 


EABERS Too” 

!' Draw the side view. 
MOVE 6,22 

IDRAW 36,0 

IDRAW 0,13 

IDRAW - 
IDRAW 0 
IMOVE 6 
IDRAW 0, 
IMOVE 2 
IDRAW 0 
! label Side 


LABEL onde 
{ Label the picture. 


MOVE 11435 

CS eae 

PACE SUCTANT Vel Saeunes 
' 


Put on the nozzle. 
i 


MOVE 14,22 
IDRAW 9,2 
IDRAW .25,0 
IDRAW 9,-2 
q 


IMOVE 0,-2 
Coie 
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NO 
& 


—~ OW GNM LS 
eee 


LABEL “nozzle” 
MOVE 14.11 
CSiZe 3 

LAseL “6° 

' 


Indicate the direction of flow. 


MOVE 65.11 

- IDRAW 2,0 
IDRAW eee = 5 
DINAN D445 
LDR A a WO 
mOwes O27 
IDRAW 2,0 
DORGINN=..o5% sO 
TDRAA Se. 5 
DRAM Tue 
SUBEND 
' 


: 
' 
! 
' 
' 


SUB Calibrate(FilenameS) 
! 

OPTION BASE 1 

Dine Coet (12) 


Calibrate the positioners on the 
milling machine movement. 

(intoesdisk, write out the calibration 
coefficients and the hard boundaries 
that must be observed! 

Has file will) be called “Metor scoet-. 


A. Position the probe volume at the 
wall of the tip of the nozzle. This 
position is (0,0,0). All readings 
will be in inches. Read all three 
potentiometers. Ask the user for the 
nozzle outer diameter and compute the 
zero position. Ask the user for the 
milling machine readings. 


B. Next, move the bed to some new posi- 
tion using the override switches. 
Take readings from the pots and ask 
for the milling machine readings. 
Compute the calibration coefficients. 


C. Move the bed to each of the extremes 
in the X, Y. and 7 directions using 
the override switches and have the 
user tell the computer when each of 
these boundaries are hit. Enter each 
of these onto the disk file, 


~~ 8 om em om Om ee a a ep fe om 8 om OR em Om Om ee om om oom som te + we 


(ik 


1220 


' 

[230 ! DD. Disk file “Motorizecoer = 

1240 ! 

1250 ! 1. x _zero, x_slope 

1260 ! 2. y_zero, y_slope 

tz70> -? 3. tere.) » sileoe 

fes0. ? 4, x_min, x_max 

12390 3! 5S: “ymin. yo max 

1300 =! Sze bile eer mice 

r3to0 

le20) BEE 

feo PRINTER IS 1 

1340 GCLEAR 

Eee GUTPUT 2 USING “#6 =255. 75 

1 ! 

1370 PRINT “I. POTENTIGMETER CALIBRATIGN: oe 
1380 PRINT “ “ 
1390 PRINT * NOTE: 1. Probe must be horizontal” 
1400 PRINT * 2. Arm must be parallel to ™ 
1410 PRINT ° the bed axis is 
1420 PRINT ” 3. ALPHA = BETA < 
1430 PRINT “ : 
1440 PRINT " A. Using the override suitches, ¥ 
1450 PRINT ” position the probe volume at . 
1460 PRINT °* the outer wall of the tip of 2 
1470 PRINT ° the nozzle. " 
P4s0 eR ae - 
1490 PRINT “ B. I will need the norzle G.0. and" 
1500 PRINT " the milling machine positon. ot 
1510 PRINT ‘ 
1520 YPRONT {Hit <eont> . > 
1S$30 PAUSE 


1540 "Beer tsc0e. | 

1850 INPUT “1. Nozzle@@8D. “inehes)? {Nozzle od 
1560 “eeeeezoc0re. 

1570 INPUT "2. X Cind. (+ eine ft lame) ee 

1580 “BSR e2s ours 

1590 INPUT “3. y Cin), (* along flume to the righee.,(-! 
1600 (BEEP cont) 

1610 INPUT "4 92 Cin) (3 custard? : 2am 

1620 ! 

1630 CALL Read_pot("X" ,Vx_1) 

1640 CALL Read_pot("Y",Vy_1) 

1650 CALL Read_pot("Z",Vz_1) 

1660 ! 

1670 ! 

1680 CALL Clear_screen 

1630 PRINT “C. Move the milling machine to a new" 
1700 PRINT “ POSttl One ihm ley ea ceteds ts > i 
1710 PRINT ” inches in each direction. a 
1720 BEEP S32 00e | 

1730 INPUT xX. y. Zim inehes? <3 «eee ee 

1740 CALL Read_pot("X" ,Vx_2) 

1750 CALL Read_pot¢("Y",Vy_2) 

1760 CA Readupot (12.7, nme 

77.0) 

_— ' C. Calculate the calibration coefficients 
1730 

1300 X_zero=*-((Nozzle.od/2)4+Vx_ I=CCX_ 2-xk 10/7(Ve2-Vxe 1D? 
1810 X_slope=(X_2-X_1)/(Vx_2-Vx_1) 


TZ 


1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1230 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2030 
2100 
2110 
2120 
2130 
2140 
2\5e 
2160 
2170 
2180 
2190 
2200 
2210 
ecau 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2oau 
230 
2340 
2390 
2360 
2370 
2380 
e300 
2400 
2410 


¢ 
Y_zeroz-Vy_1*(¥_2- 
! _slope=(Y_2- al De CU ye i) 
7 zero=-Vz_1#(7_2-Z_1 
Z_slope=(Z_2-Z_1)/(Vz_2-Vz_1) 
' 


Girection. 


Talley (ye a-Vyo |) 
(V 


Ore e471) 


' D. Find the physical boundaries for each 


CALL Clear_screen 


PRINT “1. Move the milling machine to the 


PRIN] minimum value of ‘x’ 


. Hrer<cont>. 


PRINT “(Away from the flume, backwards} 


PAUSE 
CALL Read_pot("X",V) 
X_min=X_zerot+tX_slope*V 


PRINT "2. Move the milling machine to the 


PRINT ° maximum value of ‘x’. {Towards 
PRINT the flume}. Hit <cont>. 
PAUSE 


CALL Read _pot("X",V) 
X_max=X_zerotX_slope*#V 


PRINT “3. Move the milling machine to the 


PRINT ° minimum value of ‘y’. 


Hit <cont>. 


PRINT “{To the lett along the flume} 


PAUSE 
CALL Read_pot("Y",V) 
Y_min=Y_zerot+Y_slope*V 


PRINT "4. Move the milling machine to the 


PRINT "' maximum value of ‘y’. 


PAUSE 
CALL Read_pot("Y",V) 
Y_max=Y_zerot+Y_slope*V 


Hit sscont>. 


PRINT “5S. Move the milling machine to the 


PRINT °“ minimum value of ‘2’. 


PRINT “{Downwards}" 
PAUSE 
CALL Read_pot("Z",V) 
Z_min=Z_zerot+Z_slope*V 


net, <eont>. 


PRINT "6. Move the milling machine to the 


PRINT: maximum value of ‘2. 


PAUSE 
CALL Read_pot("Z",V) 


Z_max=Z_zerot+Z_slope*V 
‘ 


Hate <cont>. 


' E. Write out the file "“Motor_coef”™ 


t 


ON ERROR GOTO Purge_file 
Reenter: CREATE BDAT F:ilenameS,1 
ASSIGN @File TO Filenames 

' 


~ QUTPUT @File:X_zero,X_slope 
QUTPUT @File;Y_zero.Y_slope 
OUTPUT #File:Z_zero,Z_slope 
' 


OUTPlerle -Xemin. x max 


OUTPUT @File:Y¥_min, Y_max 
OUTPUT @File;Z_min,Z_max 


Tes 


2421 
2430 
2440 
2450 


2490 
2500 
2510 
2520 
25a 
2540 
e5a0 
2560 
2570 
2580 
2550) 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
27 
2720 
2730 
2740 
2750 
2760 
2170 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
28390 
2300 
Zon 
2920 
2930 
2940 
23950 
2960 
23970 
2980 
2390 
3000 
3010 


! 


ASSTCNMSRiale Time 
SUBEXIT 


Purge_frle: PURGE Filenames 


GOTO Reenter 


t 


SUBEND 

' 

! 

! 

! 

! 

’ 

SUB Read_pot(DirectionS, Value) 
’ 

' Read one potentiometer and return a volt- 
! age, 

' 

' R3 -- Pot X 

' RG -- Pot Y 

’ RS -- Pot Z 

' 

' 


Directin$="X" THEN Relay=3 
Direction$S="Y" THEN Relay=4 
Direction$S="Z"" THEN Relay=5 


mi rir 
nT 


om = 


HULPUTL 7223. Deyo 
DUPPUT (723-508 71 Relay acre 
i 


OUTPUT 723;"IP.3T” 
ENTER 7240) Value 
! 


' 
SUBEND 
! 


! 
! 
! 
! 
! 
SUB Clear_screen 
Clear theo Cr ie 
’ 


OUTPUT 2 USINGR eG 3255, / > 
GCLEAR ‘ 

SUBEND 

' 


UB Motor(Direction$,RotationS) 


Turn on the motor in the requested direc- 
tion (x,y,z) with the requestion rota- 
tion (CH. Gen 


ee ey (9M ~—— o~m om o-m om « 
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3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
S110 
3120 
S120 
3140 
3150 
3160 
3170 
3180 
sro 
3200 
3210 
3220 
3230 
3240 
Sao) 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3350 
3340 
SoU 
3360 
e370 
3380 
3330 
3400 
3410 
3420 


3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
Jac 
3530 
35.40 
3550 
3560 
sore 
3580 
3530 
3600 
3610 


DirS=Directions$ 

RotS=Rotations 

IF DirdS="X" AND RotS=2"CW" THEN Lbit=2 

IF DirS="X" AND RotS2""CCN" THEN Lbit=] 
IF DirS=""Y" AND RotS="CW" THEN Lbit=4 

IF DirS="Y¥" AND Rot$="CCW"’ THEN Lbit23 
IF DirS="Z" AND RotS="CCW" THEN Lbit25 
IF DirS="Z" AND RotS2"CW"" THEN Lbit=6 

' 


! 


OUTPUT 723;"0P,7,";27*Lbit;"T" 
SUBEND 
t 


WY a O- oo am Om sp = 


SUB Motor_stop 
Stop all motors! 


OUT EUs CUR Ol. 
SUBEND 
; 


_—— s—m «—_ 


! 
: 
f 
! 
‘ 


SUB Retrieve_coef (Coef(*) ,FilenameS) 
! 
OPTION BASE | 


Retrieve the potentiometer calibration 
coefficients from a disk file called 
“Motor_coet™. Place these in an array. 


ASSIGN @File TQ Filenames 
’ 


o~—™ C= 4p OR om 


"3430 Epi ii a Sas 2 
ENTER @File;Coef(I),Coef¢I+1) 

NEXT I 

t 


ASSIGN @File TQ * 
SUBEND 
! 
' 
: 
: 
: 
: 


SUB Label_point(X,Y,2,Symbol$) 


! 
' Label a point on the Flume diagram 
' using the symbol specified. 

‘ 


PRINT TABXY(1,13).°(X.Y.Z) inches = "; 
PRINT SUSING “DOD .DDD 427.2 


7S 


t 
IF SymbolS="X" THEN 
t 


MOVE 14+Y,24+72 
IMNVE -.2,.-.2 
IDRAW 0,.4 
IDRAW .4,0 | 
IDRAW 0.-.4 | 
IDRAW -.4,0 

' 

MOVE 14+Y,11+X 
IMNVE -.2,-.2 
IDRAW 0,.4 
IDRAW .4,0 
IDRAW 0,-.4 
IDRAW -.4,0 

END IF 
t 


IF Symbol$=""+" THEN 
f 
MOVE 14+Y,11+X 


IMOVESG, .2 
IDRAH -.2,-.4 
IDRAW .4,0 
IDRAW -.2,.4 

f 

MOVE 14+Y,24+2 
GE See 
IDRAW -.2,-.4 
IDRAW .4,0 
IDRAW Src a 

END IF 


SU8END 
i 


ea ee ee ee 


SUB Move_ldv_to(X,Y,Z,Length_acm,Length_probe,Angle_arm,Angle_probe,(C< 
DEG 
' 


OUTRUN 72s Cet. 

GUIPUT 7233°S 3 372 cope 
PRINTER IS 1 

UN KEY 0 LABEL “ABORT CALL Stop_all 
OPTION BASE 1 

Tolerance=.006 


Move the probe to the position indicated 
in (inches) relative to the nozzle tip. 


Dimensions are in inches and degrees 


A.1 Load in the calibration coefficients. 


~~ 8 ew ee 


X_zero=Coef (1) 
X_slope=Coef (2) 


rae 


4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4220 
4330 
4340 
4350 
4360 
4370 
4380 
4330 


4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
44930 
4500 
4510 
4520 
4530 
4540 
4350 
4560 
4570 
4580 
4530 
4600 
4610 
4620 
4630 
4640 
4650 
4660 


Y_zero=Coef (3) 
Y_slope=Coef (4) 
Z_zero=Coef (5) 
2_slope=Coef (6) 
X_min#Coef (7) 
X_max=Coef (8) 
Y_min=Coef (9) 
Y_max=Coef (10) 
Z_min=Coef (11) 
Z_max=Coef (12) 

' 

: A.2 Move the probe 


CALL Probe_moe(Angle_probe) 


UPC oUoOesS ) (X..1 22) Lie wirthanethe per 
' mitted boundaries? 
' 


IF X>X_max OR X<X_min OR Y>Y_max OR Y<Y_min OR Z>Z_max OR Z<Z_min THEN 


Beer see S 
BEEP 2000,..5 
BEEP 1700,.5 
Beer ecu a5 
PRINT “Desired point is out of range!” 
SUBEXIT 
END IF 


C. Find out where the probe is now, draw 
the flume on the CRI, and label the 
desired position. 


1. Sound warning, movement immenent! 


ee oe ee ee ee) 


CALL Clear_screen 
PR Eerie ker 1202 rOVENENT UF MILLING MACHINE IMMINENT!!!" 
FUR I= 116-4 
BEEP IZ00,. | 
BEE aie eeet 
Beer 2200,.1 
BEEP a Ouro 
NEXT 1 
CALL Clear_screen 
' 
UTR Ue eee Cee) ‘CLEAR RELAY CARD. 
CALL Position("X",X_actual ,Valu_shaft,Length_arm,Lenqth_probe.Angle_arm. 


Angle_probe,Coef(*)) 


4670 


CALL Position ("Y",Y_actual ,Valu_shaft,Length_arm,Length_probe,Angle_arm. 


Angle_probe,Coef(#)) 


4680 


CALL Position("Z',Z_actual ,Valu_shaft.Length_arm.,Length_probe.Angle_arm. 


Angle_probe,Coef(*#)) 
. 


4630 
4700 
4710 
4720 
4730 
4740 
4750 
4760 


' 
CALL Draw_flume 


Valea be lepoint cA actual. | acttal,2 actual, +") 
Cn Gabeltooint(X,1,2, >) 


' D. Move each motor to bring the error 
between actual and desired position 


7F, 


4770 ' into tolerance. 


4780 ' 

4790 X_old=X_actual 

4800 X_node: Xerror=X-X_actual 

4810 IF ABS(Xerror)>Tolerance THEN 
4820 IF Xerror>0 THEN Rot$="CCh" 
4830 IF Xerror<0 THEN RotS="Ch" 
4840 CALL Motor("'X" ,RotS) 


4850 CALL Position("X",X actual ,Valu_shaft,Length_arm,Length_probe,Angle_as 
Angle_probe,Coef(*)) 


4860 CALL Plot_path(X_old, Y_actual ,Z_actual ,X_actual , Y_acttaweeeuee 
al) 

4870 X_old=X_actual 

4880 GOTO X_node 

4890 END IF 

4900 ! 

4910 ! 

4920 Y_old=Y_actual 

4930 Y_node: Yerror=Y-Y_actual 

4940 IF ABS(Yerror)>Tolerance THEN 
4950 IF Yerror>0 THEN RotS="CCh" 
4960 IF Yerror<0 THEN Rot$="CH" 
4970 CALL Motor €’Y jRotS> 


4980 CALL Position("Y",Y_ actual ,Valu_shaft,Length_arm,Length_probe,Angle_a 
Angle_probe,Coef(*)) 


4990 CALL Plot_path(X_actual,Y_old,Z_actual ,X_actual , Y_actual., ae 
al 

5000 Y -old=7 actual 

5010 GOTO Y_node 

5020 END IF 

S0s0: ! 

5040 ! 

5050 Z_old=Z_actual 

5060 Z_node: Zerror=Z-Z_actual 

5076 IF ABS(Zerror)>Tolerance THEN 
5080 IF Zerror>0 THEN Rot$="CW" 
9090 IF Zerror<0 THEN Rot$="CCh" 
5100 CALL Motor("Z" ,Rot$) 

Sr0 CALL Position("Z",Z_actual ,Valu_shaft,Length_arm,Length_probe.Angle_a 
Angle_probe,Coef(*#)) 

5120 CALL Plot_path(X%_actual ,Y_actua!l ,.Z_old,X_actual , Y_act@eue eee 
al) 

S10 Z_old#=Z_actual 

5140 GOTO Z_node 

Soo END IF 

iGo! 

Sioa 

5180 CALL Motor_stop : 

Soo"! 

5200 ! 

5210 FOR I=t TO 4 

5220 BEEP 2400, .2 

S250) BEEP wae Ulyere 

5240 NEXT I 

S250 1 

5260 SBEND 

S270 

5280 ! 

5290)! 

5300 ! 
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So Ole 
ase Uon! 


S330 SUB Pos:ition<(Direction$S,Value,Valu_shaft,Length_arm,Length_probe.Angle_arm., 
Angle_probe,Coef(*)) 


5340 ! 

5350 OPTION BASE | 

5o50. «DEG 

5370 

5380 “Return the position (inches) for the 


! 
' 
ooo ! appropriate direction relative to the 
' 
' 


5400 nozzle tip. 

5410 

ee CALL Read_pot(Direction$,Vol tage) 
5430 ! 

5440 X_zero=Coef (1) 

5450 X_ slope=Coef (2) 

5460 Y_zero=Coef (3) 

5470 Y slope=Coef (4) 

5480 Z_zero=Coef (5) 

5490 Z_slope=Coef (6) 

5500 ! 

$510 ! 

5oc0 DirS=Direction$ 

S500 ! 

5540 IF DirsS="X* THEN 

S300 Valu_shaft=X_zero+X_slope*Voltage 
5560 Value=Valu_shaft-Length_arm*COS(Angle_arm)*TAN(45-Angle_arm/2.0) 
S570 END IF 

5580 ! 

55350 IF Dirs="Y° THEN 

5600 Valu_shaft=Y_zerot+Y_slope*Voltage 
5610 Value*Valu_shaft+Length_probe*(1.-COS(Angle_probe))+Length_arm«COS ( 
Angle_arm) 

5620 END IF 

5630 ! 

5640 IF DirS="Z° THEN 

5650 Valu_shaft=Z_zerot+Z_slope*Voltage 
5660 Value=Valu_shaft-Length_probe*SIN(Angle_probe) 
5670 END IF 

5680 ! 

S690 SUBEND 

S700, A! 

Slane! 

Seal! 

5740 ! 

57/510 ‘ 

5760 SUB Stop_all 

5770 ! STOP ALL MOTORS AND QUIT 

S7ao ! 

S7g0 GCLEAR 

5800 CALL Motor_stop 

$810 PRINT “MOTOR CONTROL ABORTED!!!" 

5820 PRINT “CHIT <CONT> TO CONTINUE)” 

$830 PAUSE 

59840 SUBEND 

Sogou ! 

5860 ! 

S370 ! 

53880 ! 


To 


9390 
5900 
aoa 
5920 
$930 
5940 
5950 
5960 
O97 0 
5980 
Saou 
6000 
6010 
6020 
6030 
6040 
6050 
6060 


! 
! 
SUB Plotepathn(kl. 1132 leke. eee 


' 
' Plot the path of the probe on the flume 
! diagram as the motors move the bed. 
’ 
' Lower plot followed by upper plot. 
' 
PRINT TABXY(1,13),°«(X,Y,Z) inches = "; 
PRINT USING “DO. D0D -% 272722 
MOVE 14+Y1,24+Z! 
DRAW 14+Y2,24+22 
) 


MOVE 14+Y1,11+X1 
DRAW 14+Y2.11+X2 
’ 


SUBEND 
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e 
a ee ee 


¢ 
y 
’ 


' 
’ 
! 
' 
¢ 
! 
! 
' 
! 


T_SUBS 


SUB T_couple(Temperature,ChoiceS,ScaleS,No_readings,Stdev) 


SUBPROGRAM T_COUPLE 


py Brill Cwilbreth 
19 April 1984 


PURPOSE: This program is designed to 
read type T or E thermocouples and 
return the actual! temperature to 
thescalling routine. 


Temperature -~- Temperature from the 
thermocouple in deqrees F or C from 
the thermocouple identified by 


“Choiees . 
ChoiceS -- Thermocouple choice, current- 
ly ee cnnOLeiN aa NOZZLE, or, PROBE”. 
Seales) )-- Fleror Fahrenhert er “C” for 
Celsius, or= = ior histogsam in “C’. 
No_readinqs ~-- How many readings of the 
Same thermocouple should the routine 
take? 
Stdev -- The standard deviation of the 


temperature for the indicated number 
of readings in the units given by 
Scales =. 


1. Open all relays and initialize the 


A/D converter. 


WORE 3 CC. 11 
OUTPUT RRES Shoo yo lS 25512T" 
' 


OUTPUT 723;"0P,1,0T" 


2. Close the chosen relays. 


a. Ambient [ -- Type T., relays’ &,8. 
b. Nozzle T -~ Type T, relays 7,8. 
G.) Prove (= -) lvoe Ere lay. 2. 


IF ScaleS="H" THEN 
Mriet a histogram using Celsius. 
HistoqramS="YES" 
SealeS= ¢” 
Ese 
Histoqram$="NQ” 
END iF 
? 


IF ChoiceS="AMBIENT" THEN 
ioeeceeT 


eta 235 UB. 1 Gali 85 1" 
Dee 


onl 


DWOManNNAMNLSWON— 
SeDodDd0deo0deO 


IF Choice$="NOZZLE" THEN 


Types="T" 
DUTPUT 723: “OBeie7 i cee 
END IF 
f 
IF Choice$2"PROBE” THEN 
Type$="E" 
GUTPUT 723: UB, la 2 ie 
END IF 
’ 
’ 3. j(dTake an A/D conversion and convert 
! into temperature. 
! 
Sum=0 
Sum! =0 
t 
WAIT 1 


IF Histogram$="YES" THEN GOSUB Set_up_histo 
‘ 


FOR T=! TO No_readings 
OUTRUN 26 hes 3 |. 
ENTER 72301;A 

IF TypeS="T" THEN 
A=A/1000 

END IF 

IF TypeS="E"" THEN 
A=8/1000 

eel 

' PRINT "V(mV) = ":A 

GOSUB Convert_t 
BEER Am Nouriel 
Sum*Sum+A 
Sum1=Sum1+A#A 
' 


IF Histogram$="YES" THEN GOSUB Plot_point 
¢ 
NEXT I 


Temperature=Sum/No_readings 
' 
! 
IF No_readings=! THEN 
S tdev=0 
por 


: S tdev=SGR(ABS( (Sum! -No_readings*Temperature 2)/(No_readings-!))) 
ND IF 


' 
{ 4, OQpen all relays. 


! 
GUTRUTS 7237, OP leiOlle 
' 
SUBEXIT 
' 


! 
' 
Convert_t: ! 
' This subroutine converts (mV) 
' from a thermocouple into Temp- 
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1210 ! erature. 
1220 t 
124051 
[240 IF Types="T" THEN 
1250 A=2.5661297E+1«A-6. 1954869E-1#A*A+2,.2181644E-2*A~3-3.55009E-4*A~4 
1260 END IF 
(270)! 
1280 ! 
1290 IF TypeS="E" THEN 
1300 A=1.7022525E+1#A-2.2097240E-1*AA+5 . 480931 4E-3#A°3-5.7669892E-5+A~4 
loo END LF 
' 


1320 

1330 ! Fix the scale. 

1340 ! 

1350 IF Scale$=""F" THEN A=1,8#A+32 
1360 RETURN 

fers! 

iseo, 

i38¢! 

1400 ! 

1410 ! 

1420 ! 

1430 Set_up_histo: ! 

1440 ' 

1450 ' Set up a histogram of temperature 
1460 ! versus number of counts. 
1470 ' 

1480 ' 1. Zero out the Height(*) array. 
1430 ! 

1500 DIM Height(203) 

1510 ! 

[S20 FOR I=! TO 202 

1530 Height (1) =0 

1540 NEXT I 

1550 ' 

1560 { 

20 DUNE SDEVICE 1S 701 

1580 GINIT 

1590 GUTPUTs2 USTNG “er .B°:255,./75 
1600 GRAPHICS ON 

1610 FRAME 

1620 WINDOW -100,100,-10,100 

1630 MOVEs—65,,92 

1640 CSiZe a7 

1650 LABEL “TEMPERATURE HISTOGRAM" 
1660 HeeS 257,110 ,0,0),4),5,5 

1670 PEN =] 

1680 MOVE 0,-10 

1690 DRAK 0,0 

1700 MOVE 0,90 

1710 DRAW 0,100 

1720 PEN 

1730 ' 

1740 ' Take 10 temperature readings to get 
750 ! the scale. 

1760 ! 

1720 Sum=0 

1780 Sum! =0 

1790 FOR =e ag 

{S00 OUTRUN 7232 ae. T™ 
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1810 
13820 
1330 
1840 
1850 
1860 
1370 
1880 
1890 
1900 
1910 
1920 
1330 
1940 
i s0 
1960 
1970 
1980 
i220 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2150 
2170 
21380 
2120 
2200 
2210 
cac0 
2200 
2240 
220 
2260 
2e7U 
2280 
2290 
2300 
2310 
2320 
can 
2340 
Zoo0 
2360 
eo70 
2380 
Zoo0 
2400 


ENTER 7230958 


NEXT [ 


FOR I=1 TO 100 


OUTPUT 7Ageer IP 3st” 

ENTER 72301:A 

IF TypeS="T" THEN 
A=A/1000 

END IF 

IF TypeS="E° THEN 
A=A/1000 

END IF 

GOSUB Convert_t 

Sum=SumtA 

Sum1=Sum1+A#A 


NEXT I 


i 
: 


‘ 
5 


T_mean=Sum/100 
Sd=SQGR(ABS(Sum!-100*T_mean “2)/99) 
' 


f Change the window to extent from 
! -4*Sd to +4#Sd. 

' 

PVs = Nieige 

CSIZE 4 

LABEL USING "DDD.DD";T_mean 

MOVE ogee 

Oe eee 

{ 


MOVE yea sr 

LABEL USING “DDD.DD";T_mean-3*Sd 
MOVE GSee7 

LABEL USING “DDD.DD";T_mean+3*Sd 
' 


, 


MOVE 5,47 
LABEL "SO" 
¢ 


f 3. Calculate the window temper- 


! ature interval. 
' 


Intenval=4-sd) 100 
? 


T_min=T_mean-4*#Sd 
T_max=T_meant+4*#Sd 
' 


WINDOW T_min,T_max,-10,100 


Plot pOnmt. ee 
Q 


Plot each temperature paint as 
received on the histogram. 


req=((A-T_min)/(T_max-T_min))>#*3000+1000 
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2419 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
20 
2320 
2530 
2540 
2200 
2260 
2570 
2580 
2550 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2630 
2700 
Zr 
2720 
2730 
2740 


BEEP Freq,.01 


=i) 


t 


FOR T=T_min TO T_max STEP Interval 


J=J+1 

i Asad) THEN 
Height(J) Height (J) +1 
GOTG Continue 

END IF 


NEXT T 
’ 


: 
' 


Draw the point. 


Continue: ! 
‘ 


oe em FP te 


RETURN 
f 


CS) «~m ome cme ome sw 


UBEND 


MOVE A.Height(J) 
CSIZE 4 
ABEL) * = 


Label the number of points and 
the temperature. 


Pe Ue nba C273), 0CC) 2 A 
PRINT TABXY(2,5);"“Sample #";1 
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: Tea 

f 

'This program records the output from the 
'ambirent, probe and nozzle themocouples 
'for calibratren curposes. (1 nemeata cis 
‘printed as well as transmitted to a data 
i rer 

\ 


'Data 1S recorded in the following order: 
! a. Ambient temperature 

b. Probe temperature 

' c. Nozzle temperature 

' 

'Place the ambient, probe and nozzle 
'thermocouples into a bucket of warm water, 
fexecute this program, then add ice. 

' 


: 


' 

LOADSUB ALL FROM ‘'T_SUBS” 

, 

' Identify the BDAT file 

' 

INPUT “NAME OF FILE WHERE DATA IS TO BE STORED?" ,FilenamelS 
' 


Identify the number of data points desired 


' 
INPUT "NUMBER OF DATA POINTS DESIRED?” ,Nitems 


9 
' Identify the number of samples per 
' thermocouple per data point 


INPUT “NUMBER OF SAMPLES PER DATA POIN]?” Sample 


> 


' Initialize a counter 
’ 

at 

! 

! ‘Create va datay tale 


Records=(Ni tems*8*3/256) +2 
CREATE BDAT Filenamel$,Records 
ASSIGN @Filet TO FilenamelS 

' 


' Take data and print results 
’ 


CALL T_couple(T_ambient, AMBIENT”, “C’,Samole,Sd) 
BEEP Jl=So0e0> 

PRINT "“TAMB,SD=";:T_ambient,Sd 

CALE T_couple<), PROBE’ ."C “Sameve. stacey? 

BEEP T#50,.05 

PRINT “TPROBE .SD="-1.St- dev 

CALL T_couple(T_nozzle,"NOZZLE”,"C",Sample,Sd) 
Beare l*shees 


PRINT ~“1N@Zi S0= Iinezzre cu 

’ 

! Send data to the data file 

i 

OUTPUT YFilel:T_ambient,T,T_nozzle 


! Test to see if finished 
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. 


bt 


F t=Nitems THEN 
GOTQ 750 
eS 


i) 


a[l+] 
OTO 460 
EMD fF 

’ 


Gyr 


Elosel the data trle 
' 


QUTPUT 3Filet:-100 
ASSIGN @Bilel TO - 
4 


! 
‘ 
PRINT “ALL DONE” 
BEEP 
BEEP 
=ND 


Alert the user...the job is completed 
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6. PROBE CAL 


: 
! 
f 
! 
! 
! 
: 
! 
: 
‘ 
! 
: 
: 
: 
f 
f 
: 
: 
! 
: 
: 
! 
! 
: 
: 
: 
! 
: 
! 
! 
t 
! 
' 
! 
! 
f 
! 
: 
! 
: 
: 
! 
! 
: 
: 
: 
' 
! 
! 
! 
: 
! 
{ 
! 
' 
' 
: 


PRuUBescee 


This program will provide the data 
necessary to calibrate probe deflection 
if used in the following manner: 


da. 


Suitch OFF probe actuation power 
at the probe oOase. 


Suing the probe arm out of the 
tank such that the probe 1s 
accessable. 


Attach the calibration panel to the 
probe assembly, takimaqmeare Notte 
damage the glass probe. 


Run the program...it will ask for 
the desired position in mV. The 
following guidelines apply: 


tr. If it is desired @tomieuer 
the probe, type the extreme 
value 9000 (anything »>4600 
will work, but 9000 is a 
quick and easy number to 
enter). 


its Lf tt tS desired to raise 
the probe, enter the extreme 
value 390 (anything less than 
940 will work). 


._The program will next ask which bit 


tis selected to be “high". The 
following guidelines apply: 


i. If it 16s desired to lower the 
probe, enter /7. 


it. If it 1s desired to raise the 
probe, enter 8 


Switch UN probe activation power and 
when the probe reaches a desired’ 
degree of deflection, switch the 
activation power OFF and record the 
mV value printed on the screen. 


Repeat steps (d) through (f) until 
sufficient data is collected. Load 
thrs data into the programe tule, re 
and request a second order fit 

(let X=angle and Y2mV). 

Enter the coefficients derived by 
this program into the appropriate 
location in the program "PROBE_SUBS”. 


NOTE: The relays %url) be aetivateaduom 
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610 


1000 
1010 
1020 
1030 


' digital low! When the machine boots 
! up (he-9826), all relay lines are 

! high ¢+5V). The instructions below 
! will drop the voltage to zero. 

’ 


INPUT "What is the desired position? (mV)",Voltage 


Exit Desired position) (mv) = *Voltage 
' 


“!'Read the actual motor position. 


! If the desired position is BELOW the 
actual position, then tell the motor 
to move UP, 


' 

' 

' If it is ABOVE, then tell it to go DOWN. 
: | 


! 
Peo “wnles Bit DUP TUU HANT HIGH? .-bit 
' 

' 
UTP 2s, On wgene) 9! CLEAR AEE @ p70 
MUU ee co Psy 42 ebat. to 

' 


OUTPUT 723;"CC,1T" !CLEAR A/D CARD 
? 


RO 2o. Sho yces tocar Te eURMAY AND 
OU cos OR let  CLUSE ThE REEAT 
' THAT CONNECTS THE 
TAA DecAnk De Oe the 


TPOTENTLONETER 
URE 23; 1P.31° 'START A/D CONVERSION 
EhliteR e223) 7A JENTER A/D VALUE INTO A. 


OTS oA 
BEEP ABS(A),.01 
! 


! The following IF stops the motor when 
’ ait reaches it’s limits. 
' 


IF (Lbit=7 AND A>4600) OR (Lbit=8 AND A<940) THEN 
Cielo OPS 7 0 
eJelele 
eileen cee ata 

BHD IF 

GOTO 890 

END 
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7. MOTOR_CAL 


19 ' MUTOR_CAL 

20 ' 

30 ' PURPOSE: Calibrate the potentiometers 
40 ! used with the motors on the milling 
50 ! machine. 

60 ! 

70 ! 


30 OPTION BASE 1 
90 DIM Coef (12) 
f 


110 DUTRUT 7233 "SF 333 one 
20 LOADSUB ALL FROM "MTR_SUBS" 
130 LOADSUB ALL FROM “PROBE_SUBS” 
140 FileS=""motor_coef” 


150 ! 

ee CALL Calibrate teiles) 

| : 

130 CALL Retrieve_coef (Coef(*) ,Files) 

Lou BEEP 2400"3 

200 ! 

210 ! 

220 PRINT “What (x,y,z) position do you wish to” 
230 PRINT " move to? (inches relative to noz-" 


240 PRINT " zle)" 
2530 INPUT SX, 142 
f e 


260 

270 CALL Move_lIdv_to(X,.Y,2Z,Length_arm,Lenqth_probe,Angle_arm,Probe_angle,Coe 
*)) 

280 ! 

e230 BEER 

300 GOTO 280 

310) END 


20 


8. LOAD XYZ 


10 
20 
a0 
40 


POAD eK ie 

by Bill Culbreth 

30 April 1984 

PURPOSE: ‘This program will allow the 


user to enter desired (x.y,z) posi- 
tions of the milling machine relative 


to the nozzle. The values will be 
stored on disk to be utilized later 
by MAIN_T. 


Et A000) 51500), 71500) 


le) oimput the file name. 


_—~- o—w ome [ J ee ee ee ee ee ee ee 


GOSUB Clear_screen 


Pron ener s | 

PRINT “Input the file name.” 

PRINT "“{ I suqgest “RUNXX’ where °XX’*" 
PRINT " aus the run number. 4 
PRINT 

INPUT Filenames 

' 


2. Begin tnputting data. 


om ome oP ome 


GOSUB Clear_screen 
' 


PRINT “Do you wish to append a previous data file?” 
INPUT Answers 
' 


IF AnswerS="YES" THEN 
INPUT “Previous file name?” ,Old_files 
ASSIGN @Filel TO Old_filed 
. 


1=0 
Koop l: ! 

T=[+] 
SNeRaeel belextt Gl). 204) 
ie Oleic = le KOT). (C1). 2 CED 
IF X¢€I)<>-100 THEN GOTO Loop!i 
ASSIGN @Filei TO * 
PURGE Old_filed 
Count=1-1 

BIESE 
Count20 

END IF 


GOSUB Clear_screen 
BEE 


t 


PRINT Sino tethemdesmrad Desitron in 
ERENT Inlenes aS. sv aoe nidel. 

PRINT 

PRullie 22 lecminate Inewt by entering -100°” 


ont 


ftoRe-mem woe oars wer ee 
COWOODNMANUNLSWN—- © 
OooQoooooOoO0O OC © 


PRINT “2. Terminate ineutwoy entering 1 00a. 
PRit for xX, Y¥oueandec 
PRINT "3. If you wish to enter an ortentation angle,’ 
PRING entier “-S999orremtacron, 0 . 
PRINT 
t 
Begin: ! 


Count =Count+] 

PRINT “Item #";Count 

INPUT “€X,Y,2) in inehes 2 4 .cCount) ,¥(Coune), 7 CCount, 
OUTPUT 7013 "1, xt Z="sCoume, xX CCount), ¥( Geant) {2 Coin 
IF X¢€Count)<>-100 THEN GOTO Begin 


End_count=Count-1 
' 
All data points have been entered. 


{ 

! 

! a. Set up new softkeys. 
f b. Explain softkeys. 


GOSUB Clear_screen 


PRINT “SOFTKEY LABELS:” 
PRINT “ Q == WRITE data to diske] 
PRINT 2 -=-"EDP) out bad dataw 


PRINT * 4 -- HARD copy the data on printer.” 
PRINT “ 6 == LIST data on the Chae 

PRINT " 8 ~-~ STOP terminates the program.” 

! 

ON KEY LABEL “WRITE GOSUB Write_data 


0 

2 LABEL “EDI GOSUB Editmdata 
ON KEY 4 LABEL “'HARD’ GOSUB Hard_copy 

6 LABEL “LIST” GOSUB Listidata 

; LABEL “STOPS GOTO Terminate 


Clear_screen: i 
t 
! Clear the CRT display. 
' 
OUTPUT 2 USING w#,BN5 255.75 


GCLEAR 
RETURN 
! 


<< =m om om oc 


! 
Hard copy: ! 
' 


Print out all data to thesemunter. 


PRINTER IS 701 


on 


1240 
1220 
1220 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1o30 
1340 
1350 
1260 
1370 
1380 
1290 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
SZ 0 
1330 
1540 
1550 
1560 
7 0 
1380 
1530 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1630 
1700 
here 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1730 
1800 


PRINT " Count X€inches) Y¥<count) 
PRINT " 
PRINT 
POR [=i 10) End count 
PR EN ieee Cl) Ch 7 6.) 
NEXT [I 
? 


PRINTER IS 1 
RETURN 


' 
' 
' 
’ 
' 
List_data: ! 


' List data to the CRT. 
9 


GOSUB Clear_screen 
’ 


Pai There ameh +End count; data points.” 
PRINT 

Lari nnvehepoint do l start with?” ,Start 
INPUT “which point do I end with?”,Ena_data 
’ 


FOR I=Start TO End_data 
aN olka t eee, KGL (2) ZC] ) 
NEXT I 
RETURN 


: 

: 

! 

t 

Edit_data: ! 
’ 


' Edit data. 
9 


GOSUB Clear_screen 
' 


INPUT “Which data point do you wish to alter?” 
’ 


PRINT 

Cutie “-oOreooumnees «0% 4 (X,%.2) where:" 
Pie ie Ar tenet Fp sch ote 7 eT) 

' 


INPUT “Tyger in the mew values: ,X(I),YCI),Z¢1) 
RETURN 
' 


’ 
: 
! 
' 
Write_data: : 


UEWhlcte sdatacout’ to a file on disk. 


a 


Z7(count?)” 


I 


1810 ' 


[820 XCEnd count *1)=— 100 

1830 YCEnd counttlo== 100 

1840 ZCEnd_ count+!1)=-100 

1850 ! 

1860 Max _data=3*(End_count+1) 

1870 : 

1880 File size=INT(Max_data*8/256) +1 

1890 ' 

1900 CREATE BDAT Filename$,File_size 
1910 ASSIGN @File TO Filenames 

1920 ' 

1930 FOR T=#1 TO End_count+l 

1940 OUTPURSSRi le :X¢( 1), 1Ci)2 cD) 
1950 NEXT I 

1960 ' 

1970 GOSUB Clear_screen 

1980 BEEP 2400..3 

1380 PRINT "File “;Filename$;"" has been stored!” 

2000 ASSIGN @File TO * 

2010 RETURN 

2020 =! 

2030 ! 

2040 ! 

2050 ! 

2060 ! 

2070 ! 

2080 Terminate: ! 

2090 GOSUB Clear_screen 

A a PRINT “NORMAL TERMINATION OF PROGRAM!” 
eo 


oe 


or 


SEND DATA 
10 : SEND_DATA 
20 ' 
30 } To VAX SEBM. TRS-80. 
40 ! 
50 ' HP-9826 TERMINAL PROGRAM 
60 ' (REGUIRES BINARY ENHANCEMENT PROGRAM 
70 ! Begs id 
80 ! 
90 JUNE 30, 1982 
100 ' updated 1/5/83 
110 ! updated 1/16/84 
120 ’ 
130 Pare CULBRETH 
140 ! 
150 ' 
160 Sc29 eR So - 2 3eoms Sele et COME 9. 
78 Potirer 1S 1 'SRRINTER [S Cr. 
180 Prat Plepnuiel PRINTER 1S CRT 
190 Printer_choice=701 ! MY PRINTER IS 701. 
200 Bits=7 ' BITS PER CHARACTER 
210 Duplex=0 PULL DUPEEX 
220 Baud=300 ' BAUD RATE 
230 Computer=1 ’ ASSUME IBM COMPUTER 
240 ' 


250 OUTPUT Pr;"{(300 BAUD, IBM assumed.” 
260 CUTE Ue, seoadmene Dinarysorogram BEB first” 
2/0 QUTPUT Pr;" unless you have BASIC 2.0” 
280 GUIPUT Pr. Ser. MODEM GON <FULE DUPLEX> +" 
2390 OUTPUT Pr;:" " 
300 ' 
310 DIM Name$(200] ,Hp_f1le$(30] ,AaC1500) ,NumbS{ 30] 
320 INTEGER [send 
' 


340 CONTROL Se ,3) Saud 
Jo) CONTROL Sc,47Bits-5*4 ! BITS/CHAR & <STOP BITS. 
360 ! 

370 ' 

380 To_disk =0 
230 Datadump=0 
400 I_data=! 


410 ie 
420 2 
430 aT 
440 Et 
450 : 


460 ON ERROR GOTO Errors 
470 ON KEY 0 LABEL “Line Mode" GOTO Line_mode 
480 ON KEY S LABEL “Terminal” GOTO Terminal 
490 Oe Ke! 6 LABEL “TosGe= GOTH Pr.crt 
500 Tir eieat sabe lo Pre GUID Pripret 
aI. ON KEY 8 LABEL “DATA GOTO Data_dump 

’ 


530 ! 
540 Line_mode: ' 
550 QUimUerr: VeINeere Cem On, MODE,” 


200 Begin: STATUS Sc,10;Y ! CHECK FOR FULL BUFFER 


580 | Ieee iC yO) =0 SHENECOTO' Sein 
5600 ' RECEIVE ROUTINE. 


o> 


OwmonnwiIlS WN 
QOooooooqoooeoo 


' 
Bee ey STATUS Sc,6;A 
=A 
OUTPUT Pr USING “#,A";:CHRS(B) 
IF B=63 AND Datadump21 THEN GOTO Data_dump 


IF B=13 AND Computer23 THEN OUTPUT Pr:CHRS(13) 
GOTO Begin 


' TRANSMIT ROUTINE. 


IF Duplex=0 THEN 
IF NUM(Key$)<>255 THEN OUTPUT Pr USING “4,A";KeyS 
IF NUM(Key$) 2255 THEN OUTPUT Pr;” " 

END IF 

IF Computer#1 AND NUM(Key$)=8 THEN Key$=CHRS (64) 

' 


- 


the previous line gives an & 
' for a backspace for the IBM. 

! 

IF Computer=S AND NUM(KeyS)#8 THEN KeyS#=CHRS(127) 
' 

' 

' 


THE VAX/VMS REQUIRES A DELETE 
! SYMBOL FOR A BACKSPACE. 
' 


IF NUM(Key$)=255 THEN KeyS=CHRS(13) 


DUTPUT Sc USING "“4,A";Key$ 
GOTO Begin 


DATA FILE OUT TO THE HOST COMPUTER. 


om comm 6 fm 8m +m 


et 

Data_dump: ' : 
IF I_data=! THEN GOSUB Open_file 
} 


IF Datadump=0 THEN GOTO Begin 
IF Computer=1 THEN WAIT .3 
' wait for the slow IBM. 
BEEP 1000+RND*1500,.05 
OUTPUT “Er s AG. I data w)=-: 
DUTPUT Pr;:Aa(I_data) 
GOSUB Send_number 
IF Aa(I_data)=-200 THEN 
I_data=! 
Datadump=0 
END IF 
I_data*I_datatl 
GOTO Begin 
‘ 


! ERROR HANDLING SUBROUTINE 


Groner OFF ERROR 
gee _file=#-200 


FIRST, END OF FILE ERROR. 
iF ERRN=59 THEN 


Aa(I)=-200 


96 


1210 GOTO 2000 ! RETURN AFTER ERROR. 


1220 END IF 

| 2280) ' 

1240 IF ERRN<>S9 THEN QUTPUT Pr:"<error #'°;ERRN;" generated.>" 
1250 IF ERRN=S4 THEN QUTPUT Pr;“(FILE <“;sHp_fileS;"> ALREADY THERE:"’ 
1260 IF ERRN=S4 THEN GOTO Created 

1270 LESERRN=oo TREN TOUIPeIUT Pr. <FILE (“sHp fileS;* IS NOT GN DISK.>” 
1280 ASSIGN @File TO «* 
1290 GOTO Line_mode 

h30Ont 

130 4 

feZOee aoe) TO-CRI, 

i330 7! 

1sa0 Preocrt: Pr-{ 

iga0 GOTO Line_mode 

1360 ! 

1370) 

1380 !' QUTPUT TO PRINTER. 

Vo30 

1400 Priprt: Pr=Printer_choice 

1410 GOTO Line_mode 

1420 ! 

1430 ! 

vee ! CHANGE THE TERMINAL CHARACTERISTICS. 

1450 ! 

1460 Terminal: ! 

1470 QUTPUT Pr;" 1. Baud Rate #";Baud 
1480 GUTPUT sei Zee Bits/Char 4 Bits 
1490 GUIPUT Pr” 3. Duplex =""-Duplex 
1500 GUTPUT Pre: C1=full,O=half J" 

1510 OUTPUT Pr:" 4. Computer =";Computer 
1$20 GUTPUT Pr;" CIBM=1, VAX/UNIX=2, " 
1530 OUTPUT Pr;"™ TRS-8023, Cyber=#4, vax/vms25]" 
1S40 GUC ee = 

1550 INPUT “Change which one?” ,Which 

1$60 IF Which=1 THEN INPUT "To?" ,Baud 

1870 IF Which=2 THEN INPUT “To?",Bits 

1580 IF Which=3 THEN INPUT “To?” ,Duplex 

i5o0 IF Which=4 THEN INPUT "To?" ,Computer 

1600 IF Computer=1 THEN Duplex=0 

1610 IF Computer#3 THEN Duplex=0 

1620 IF Computer=3 THEN Bits=8 

1630 IF Computer=5 THEN Duplex=! 

1640 GOTO Line_mode 

1650 ! 

1660 ! 

1670 ! | 

1680 Open_file: ! ‘ 

1690 ' Open a file to read data from 

1700 i disk. 

ee0 ! 

1720 Datadump= 1 

1730 ! 

1740 Cie lcethisceelvedatar viwzresS: .LdavsS 

1750 IF Ldv$S=""1"" THEN 

1760 INPUT “Experiment #?" ,ExperimentS 

i770 BIESE 

1780 OUTPUT Pr:"“Data file out of HP to host.” 
1790 Nevo ples name jhowtt les 

1800 ENDeLE 


ey 


1810 ' 

1820 IF LdvS="1" THEN 
1830 He_fileS2ExperimentsSé" RESULT” 
1840 END te 


1860 ! Read the file of f of disk. 
1870 


1880 ASSIGN 9File TO Hp_files 
T=1 


1900 Check =0 

12:00 BEEP 

1920 BEER 

1930 OUTPUT Prs"{Working on file <“sHp_fileS;3">.}" 
' 


1960 ENTER @File;AacI) 
Taiz Check =Aa( I) 
1980 [=[+] 


! 
2010 ASSIGN @File TO » 
2020 Datadump=1 
ae RETURN 
' 


Z0s0) 
2060 ! 
2070 Send_number: ! 
2080 ! SEND A NUMBER UNE CHARACTER AT 
aor ! A TIME TG SHE HOS! COMPUTER. 
1 ! 
ain Numb$=VAL$S(Aa(CI_data)) 
2120 Length-LEN(Numb$) 


2140 rf ((LdvS="1") AND (I_data>!3)) THEN 
2150 Posit=POS(NumbsS,".") 
2160 IF (Posit<>0) THEN Length=Posit+2 
2170 END IF 

, 


2190 FOR I=! TO Length 
2200 Numeric2NUM(NumbS({I,IJ) 
Zero DUTPUT Sc USING "4,A":NumbS(I,I] 
2220 NEXT I 
! 


2240 QUTPUT Se USING ‘#,A";CHRS(13) 
2250 RETURN 
0 ’ 


2270 ! 


2280 ! : 
2290 END 


98 


ike 


TCAL 


ANNANNANAN 


20 


39 
50 


APEENDIX Cc 


MAINFRAME PROGRAMS 


eA 


PURPOSE: THIS PR8CGRAM PLACCS THERMOCCUFLE CALIBRATION DATA 
RECEIVED FRO™ TRE HP=—9826 MICROCOMPUTER INTO A 
MORE «ORK ACLE FORMAT FOR USE IN THE PROGRAM TFIT. 


DI“ENSTON TPROBE(200}.Tame( 200) .TNOZZ{(200} 
{=1 


READ(07.,.2) FAMA(T ) 

IF{ FamA( 1} .EO.-200.0) GO TO 20 

READ(O7.+} TFPROBE (1) 

READ( 07.5) TNUZZ( 1} 

NITEY*S=I 

f= el 

GO TO 190 

CUNTINUE 

@RITE( 5.40) 

00 30 {[F!1.eNITEMS 
@RITEfo.SOIL-TAMH( 1} -TFRPOKBE( IT} -TNOZZ(1} 

CONTINUE 

FORMAL (GX. * AMBIENT (C}*,.2X2*°PROKBE (CPL 2X.*NOZZLE ([C}*} 

FOMMATOCIixXe SeOXe2(FG6.320X)) 

STAR 

ENDO 


ae 


Z| Dean 


Let 


PURPOSE: THIS PROGR4™ PERFOR™“S A FIRST UROER CURVE FIT BY THE 
LEAST SQUARES YMETHOO FOR Tw CCLUMNS OF DATA. THE 
FIRST COLUMN LISTS VALVES OF KX AND THE SECOND LISTS 
VALUES OF F(X} OR YY. FOR THER“OQCQUPLE CALIARATIUN. 
LES x = THe AMBIENT FE*PERATLRE. AND LET Y = EL THER 
THe PROBE OR NOZZLE FEMPERATLURE aS DESIRED. 


NNNANNNNNAN 


OIMENSTON x(1S50).Yf150}.xSO(150).-xY(1E0)1.-YEST(150) 
SUMx>s0.0 

SUMxSQO=0.0 

SuUMxY=0.0 


SUMKXY=ZSUMKY 6 
SUMY=SUNYe¥Y ( 
10 CONTINUE 
MET TEMS=FLOATINI TEMS } 
A= (1 SUMYSSUMXSC-SUMXSSUMXY ISIE X LTEMSSSUMXSO-SUMxS 
Bs ( XE TEMG SSUMKY~-SUMXSZS5SUNMY JS (IX ETEMSSSUNKGO-SUMNKS 
YRARSSUMYSHILTEMS 
O0 20 TtteNtiTeEe*s 
YEST(tTpseusx( thea 
WRITE (G.2) veSsrTrtt). vlt) 
SUMMUMZSUMNUNM? (YEST (LT I-YEAR 
SUMDENZSUMDENe LY( TI) -7YA AR) 2s 
290 cCUrnerFIiINue 
RS O=SUMINUMS SUMOE®I 
aw | TE { Ae 2%) & 23 eRGO 
STOP 
END 


) 
( 
t 
SUMKSO=SUMXSU 
x 
[ 


=e 
=o 


S% 2 


“- 


100 


on 


ay 


THIS PROGRAM FITS JET FRAJECTCRY DATA TQ THE EQUATION 


THE LEAST SQUARES “EFHOO ANO COMPUTES A 


CORROLATION COEFFICIENT (REF2). DATA POINTS ARE READ 
IN FREE FORMAT FROM UNIT "7" CISK, 


JETCURV 

Cc 

C JETCURV 
C 

C 

Cc PURPOSE: 
C Z=Aasys(uey) 
e 

s 

c 

c 


10 


29 


30 
40 


DIMENSION Y(S0)-Z(50)-YY(S50)-Z2Z(50) -YSQ(SO}-YZ(S50)+ZESTISO) 


SUMY=20.0 
SUMYSYQ=0 eV 


SUMNUM=Q0 


line meme «@ ZOO 


~~ il tt<~ nm —e e 
Kee ate 


SUMYZ=SUMYZ?*Y 
SUMZ=>SUMZ+ ZZ ( 
CONTINUE 


KLTCMSZPLUATINUE TEMS } 


AAz=(SUMZSEQOUMYSQ-SULUMYSSUMYZ)S l 
BH=(XLTE“MSSSUMYZ-<SUMYSSUMZIS (XL 


A=t.0/7 4A 


Brash 
ZBHARADSUMZSXITEMS 
00 20 =teN( TES 


ZEST(t)=dB2vy ( 
wR{Trlo.e. 30)2Z( 
SUMNUM=SUMNUM & 
SUMNEN=SUMDEN? 
CONTINUE 
RSIA=SUMNUMZS SUMNEN 
ARI TE (G6 e640 Ae e RSC 


FOWMAF(EX.*2 (ACTUAL) 


XLTEMS SSUYYSO—-SUMY2F2 } 
TEMSSSULUMYSO-SUMYS=2 } 


{) 
-~Z7yARJ SH 
BAR | s22 


Sao = 


=e ti2ZeheJXa* 2 (EST =*,FIL2.8) 
FORMAT (IX. * az fseF i 2ehs 2X, *S= "eFL2e.Ge 2X 2.1 REVZ= aie aes | 
Sine 
END 


Poa 


A 


Q) 
ANANAANAANNANAANANNAANN S 


-—“AND 
Oo o 


mM SJANNNA DF N 


ANN 


GRAS 


PURPOSE: DATA TRANSFER FROM THE RP-~9EZEMICROCOMPUTOR 
TO fwHe [8™,. 


OY BILL CULBRETH 
FOR ME2410+ FALL CUARTER. 1982 


FILEOEF OS TERMINAL 
FILEODEF O6 TFRMINAL 
FILEDEF O7 OILSK “NMYDATA DATA (PERM) 


GLOKAL TXITLIA FORTMOO2 MOC2ZEER 


TYPE (NN THE ARCVE «© LINES TO MAKE THIS 
FORTRAN PROGRA™ RUN. 


OIMENSTON OATAI30N0N0O} 

{=I 

wRITFE (6.80) 

FORMAT (2X%.°*GBEGIN INPUTTING DNATA FROM TRE HOP-98B26'" 


CONTINUE 
REAUIS.+) DATA(T) 
x=[ef 
[F(OATACI-~-f£).NE.-200) GOTG 190 


NITEMS = [1 
FORMAT(2x.15.* NATA POINTS «ERE ENTERECO.,'*) 
WRIETEL6.641 NITENMS 


NOW THAT ALL OATA HAS RPEEA ENTEREO.e wREITF [T QUT ON 
O1SK. 


FO2MAT( 2K. *OATA(*® -£5e%) = * - tFl Se 5? 
t=t 

ARI TEL 7.) OATA(T I 

[=tet 

[FtTOaTat(t—-t}eNFE.~2-20C!) GOFTC 20 


ALL DATA HAS MEEN wRITTEN ONIO OISK. 


STaPr 
ENO 
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Sa LDA 


TOATa 


PURPCSE: TH(S PRCGRAM FLACES FRE SUCYANT JET DATA RECEIVED 
FrRO™ IME RmP=-G52€& MICRCCCYEUTE INTC A URE ORDERLY 
FGR WAT. {TT 4uSC CCAVERTS VAITS CF LENGTHY FRGM SI 
TO METRIC. 


ANADADNAADDAN 


OIMENSION X1[1200).-Y(2901.-ZI1Z00I1.TPROQHEIZ00[ .TaAamB8l(200}) .TNOZZIZ900I[. 
ISTOEV[200) .4™@™m[ 209) .Y*"4[ 200) -Z™™=( 200) 
(=t 
10 SEaAOl(QO7.=) xt) 
[FI xc£T1) .EC.=-200.3) GG TQ 20 
meet {)=25.4SK[1) 
READ(U7.S) vei) 
ye@([)=25.a4-5”%Y( 
READIO7.,.= Zit 
Ze@([} S25 .-457(() 
R=AO(O7.+=) TeESNE 
RFEACD(C7.S) Tansy[ 
Zs 
J 


oS) = wr om 


RFEADIO7.=) FNOZ 
REACIQ7.=) STCE’ 
N(TE™S= ( 
t= +) 
GO to 10 

29 CONT (SUE 
ex [TE(5.49Q) 
90 30 ([=teNITEVS 

eX[TEeloe-SCil -xeVMild)eyv¥vl f)-Z2Z44MI()e-TERQSEL LD. TasMGif).IncZZ(t). 

CStDOev rai} 

39 CONT (NUE 


QO F@HMSstiaskt.*k [yt , ax, ty (ye. ax. 8 2 Peet we Xt PPPOE «ICI ke 
tome lcawt (Cp ° -PKe'NCZLLe (C)* 1k. *STO Ce toe 

SO “Gea Pit xk E30 JX eATR 7 cece SK Per Fed eK eI FH 30 4X) 1 
sree 
ENO 
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6. CONTOUR4 


CONTOURS 


PURPOSE: THIS PROGRAM [5 DESIGNED TC CISPLAY BUDY4ONT JET 
TEMPERATURE DATA USING A CONTOUR PLOTTING PACKAGE 
AVAILABLE wiITH OISSPL4. TrE FOLLORING Raw OATA [IS 
RF AD FAHCM A OISK FILE? PROEF POSITIONS IN xy¥2 
COORUINATES RELATIVE TO THE ADZZLE TEP. AMbIENT 
TEMPERATURE. PROBE TEMPFRATURF A4ND NOZZLE TF MPERATURE. 
THE PROPE 4NO NOZZLE TEMPERATURES ARE TRANSFIJRMELD 
BY C4LIRARATION COEFFICIENTS 4&ND NORMALE ZEO wRT Anm@BIENT 
TEMPERATURE. THE KYZ COURCINATES ARE CONVERTED TO 
XSew COORDINATES RFLATIVE TC THe INTERSECTION OF THe 
D4T&A PLARE «a@ETH THE CEATERLIAE TRAJECTORY UF THE JET 
(THE S-AxlS £5 TANGENT fO Rte fRAYECTORY). THIS 
PROGRAM ALSO CCMPUTES THt FATE OF HEAT TRANSFER EFRON 
THe YET TO ThE AMALENT, 


NOTE: TRE FOLLOWING VALUES MUST AE INSERTED AS THE 
FIRST LINE OF CATA EIN FREE FCR™AT?32 THE TOTAL NUMBER 

OF DATA POINTS ENITEMS). TRE ACUTE ANGLE RETWEEN 

THE OATA PLARE AND HORIZONTAL (THETA) AND THE VERTICAL 
OISTANCE PETac&N TRE Z-4xX1S AND THE [INTERSECTION 

POINT CISCUSSEC ABOVE (24). THE CENTERLINE VELOCITY 

IN “74S (VEL) ANDO JET wWIlOTR [XR 4M (WIOTH), 





ANNAN ANNANANDNANANANANANNANANNH 


OIMENSION X(1O00)-VYITEOOIeTELOO) -TP(LOOIV- TAL LIOOI*TNI100).TMAT( 10.10) 
COMMON wORK( 160000) 


Cc 
Cc 
C+ = = = READ DATA FROM UNIT "7" OTSK = ¢ = 3 
Cc 
Cc 
REAO( 7.) NITEMS.,TRETA.,ZA-VELewICTH 
Cc 
00 20 [FI NITES 
READ( 7,5) SeX(CPAeVOTI°~TPIC UC) TAC TITNC TILA 
Y( ET) = CYCCVRZAIDSSENE THETA) 
c 
Cc CALIBRATION CNEFFICIENTS 
TP(TI=.-9AGBODETSFATO( LD &1. 74079665 
e TNELT PRE 6-9ISLZ2Z 34 26STN(LP*l. as7540 G6 
Tit) = (TPOCED — TALCTIISETNO TI) = TAltdd 
20 CONTINUE 
Cc 
c 30 FURMAT( 2K. °Le Xe Ve Fa*el5.3JFIS52-5) 
Cc 
C ; 
Gs = = $= DI¥¢ENSTION OF THAT = S|] |= s 
G 
Cc 
IxOtI™ = [£0 
IYO iM = LO 
G 
Cc 
Cc 
a ¥ oo: = FIEND THE MAXIMUM AND “MINEI“UM TEMPFRATURES = = & S$ 
ec 
TMIN = 100.0 
TMAX = QO.Q 
XMAK=N.Q 
YMAX=0.0 
XMIN=1L0020 
YMIN=100-.0 
Cc 
00 40 T=1lLeNETES 
IFC Tt) ~Gre TMax) TMAx = TLE) 
JFETCT) ecLTe FMENY rein s Treti 
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PeGl .XMAXTXMAXEX (11) 
LPeGTeYMAXDPYMAXSY (1) 
EPeLF -XMERIKMINEX (1) 
(Pe-LTCLYMINIDYMINEY (1) 
40 CON 


4} 
tt 
4 
4 


OE TERMINE THE MAXIMUM VALUES OF XK AND Y wITHIN A 
QUADRANT, ASSUMING THE STREAMW(SE AXIS [1S CENTERED 
IN &A SUUARE MATRIX PLANE 


AAAANM 


Ox=( ( KXMAX—XMINISFLOAT(L(IXDOIMIIZLOOO. 
DY=((CYMAX=-YM(NISFLOAT(IYONIM) 1/1000. 
QUADXS=IFLUATI( IXUIMI42.0)2CX 
QUADYSIFLOATI(IYOIM)I42.0)13CY 


tt 
+? 


> = FIND THE [NCREMENTAL AREA FOR RFEAT TRANSFER CALCULATIONS 


tt 
4? 


AREA=OXK SN Y 


= ¢ ESTABLISH THE 


4s 
4 
4 


INCREMENT SIZE FOR CONTOUR PLOTS = = = 3 
TINCR=( TMAX=TMINIZ5 225 


= = GENERATE THE GRIC = = = & 


4+ 
4+ 


CALL 
Call 
CALL 
Cace 
CALL 
CALL 


NM ONANNNA ANA ANN 


CALL 
CALL 
CAtL 
Cau 
CALL 


AOAN 


CALL 
CALL 
Cale 
CALL 
CALe 
CALL 


AM 
te 
4 


CALL 
CALL 


S = GENERATE THE 


COMPRS 

TEK614 

PAGE 8e-.-9-) 
BLUWwUP (0.75) 
PHYSORI1.2S5-1.) 
AREA2N15.5.5.-5) 


RETGHT(90.100) 
CARTUG 
INTAXKS 
XNAME (*X 
YNAME ( @wW 
C4Llu XTIC 
CALL YTIC 
CALL YAXA 


2KKR=-— 


Sw SS 


READIN( P TEMOERATLRE CONTOURS IN 


4 EVUYANT JETE*®?.-{00.224.4) 
HFADINI* «if TR A CROSSFLOWING AMHBIENTSE® -100-e262%) 


HEADINI *( 25% FLCa ARATEPS* LOC 1 2 oe) 


HEADING * (PLANE AjPE*-LIDe Le o% } 


GRAF (—-30.0.205.-0e-50.90.--30.20.052-0-230.0) 


SCOUMUNILOO0ODN ) 
ZBASEITMIN) 


CALL EPGNMATIIXOIM,IYOI™)} 
CALL GETMATIXK.Y.TeNIETES,.O) 
CALL FENOMATIIMAT,.O} 


INTERPOLATED TEMPERATURE 


MATR {xX 


SOLVE THE RATE OF HEAT TRANSFER &$ 3 S & 


AOAN 
t+ 
¢} 
¢? 
tt 


asu”=0.0 
bvO 70 [=letxXOI™ 
00 650 J2ztefYOI™ 


RNA=I1teVANDSAI3I © JY +Twabrfit.J) 


3 
“A -.IDA SHOOTE-CezIvar;,# ( 
© S$-41 290029 F-RP2 Tas 
“9%. 74007EFHILST 4aAT | 
RRKRIQUADKX-FLOATI( I-12} 
YY=QUAUYHFLOAF( (=-1) 
HADIUS=SULRIT(xAKS52 © 


OD 


>1000. 


tt 


‘arene 


4? 


60 
70 


80 


4 
4 


V=VELFEXP(—-RAOLUSSS2/wlOTHss 2) 
Q=ROSANEASVS -445S5(TMAT( IL JID E73-16) 
OSUM=0SUM+Q 

CONTINUE 


CONTINUE 
WRITE l6.40INSUM 
FORMAL (IX. *9O= §.F20.5) 


= PLUT THE CONTOURS = = = 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CAcEe 
STOP 
ENO 


CONMAK(TMAT. [XOIMeILYOIM.e TINCR) 
CONLIN(O0.*SOLICS*® -*LAAELS*.~L 10) 
CONANG(90.VN0) 

RASPLN(0.25) 
CONFUR( 1. *"LAGEFLS*.*ORAWw® ) 

HE TGA (9.05) 

RLMESS(*ExP 1C $*,.100.-594-59) 
ooT ° 

GRIDNCL.1) 
tNOPL(0O) 
DONE PL 
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APPENDIX D 
TABULATED DATA 
TAB GE wl 


ROTOMETER CALIBRATION 


Flow ml/s Std Dev (ml/s) 
LO 5.94 Oreo 2 
IES) S02 to 
20 OO 0.41 
25 a OS 0.24 
30 14.03 rss 
pe 1529 8 0.58 
40 17.41 mol 
45 19.09 0.34 
50) 20 75 Oi 3S 
55 22S O23 
60 Dee 27 0222 
65 Om Oe 20 
70 27.64 Oe 
i eos 0.24 


Lig y 


TABLE 2 


TEST RESUS 


Crossflow yvVereeiE,- -13 m/s 
Nozzle flow rate: 11.85 mis 2s 
Nozzle inside diameter: 7.144 mm 
Nozzle discharge velocity (mean): 29.558 mms 
Nozzle temperature (mean): ae. © 

Ambient temperature (mean): 29 OC 

Proude Number: 14.8 


Michaelis-Menter equation: 


Coefficient A: 264 284375 
Coerittelenrer. 1.03698254 
Plane Y (mm) od (degrees) 0 (W) 
A ee 46 2 © 
B 2 0 58 Ge 2 
c So 708 70 15.349 
D 62.889 80 18.366 
E OF ones 86 2 eS 
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